. Universidad de Antioquia 

! Medellm (Colombia) 

O 

> 
O 

:z; 

(N . 

: Study of core collapse neutrino signals 

and constraints on neutrino masses from 
a future Galactic Supernova 



in 
o 

i-^l A dissertation submitted in partial satisfaction 

I ' of the requirements for the degree 



Doctor in Physics 



by 



5^ ; Jorge I. Zuluaga 



Thesis director 
Prof. Enrico Nardi 



January 2005 



Abstract 



We study the sensitivity to neutrino masses of a Galactic supernova neutrino signal 
as could be measured with the detectors presently in operation and with future large 
volume water Cerencov and scintillator detectors. 

The analysis uses the full statistics of neutrino events. The method proposed uses the 
principles of Bayesian inference reasoning and has shown a remarkably independence 
of astrophysical assumptions. 

We show that, after accounting for the uncertainties in the detailed astrophysical 
description of the neutrino signal and taking into account the effects of neutrino os- 
cillations in the supernova mantle, detectors presently in operation can have enough 
sensitivity to reveal a neutrino mass (or to set upper limits) at the level of 1 eV. 
This is sensibly better than present results from tritium /3-decay experiments, com- 
petitive with the most conservative limits from neutrinoless double /3-decay and less 
precise but remarkably less dependent from prior assumptions than cosmological 
measurements. 

Future megaton water Cerencov detectors and large volume scintillator detectors 
will allow for about a factor of two improvement in the sensitivity; however, they 
will not be competitive with the next generation of tritium /3-decay and neutrinoless 
double /3-decay experiments. 

Using the codes developed to perform the generation of synthetic supernova signals 
and their analysis we created a computer package, SUNG (SUpernova Neutrino Gen- 
eration tool, http: / / urania.udea.edu. CO / sungweb ) , aimed to offer a general purpose 



solution to perform calculations in supernova neutrino studies. 
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Nothing is built on stone; 
all is built on sand, 
but we must build as if the sand were stone. 

Jorge Luis Borges 
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CHAPTER 
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Introduction 



Supernovae (SN) and neutrinos form one of the most interesting (and strange) cou- 
ples in Astrophysics and Fundamental Physics. The first ones are the most energetic 
explosions in the Universe after the Big-Bang while the second ones are the lightest 
and most elusive (massive) elementary particles. 

This strange relation was probably first recognized in the early 1940s P and arises 
from the fact that neutrinos are produced in huge quantities by the nuclear processes 
which dominate the evolution of a dying massive star. 

When the iron core of a massive star M > 10 Mq reaches a critical mass Mc — 
1.0 Mq (Ye/0.4)^ j2] by accumulation of the residual material produced in the "melt- 
ing" of silicon, it becomes unstable and starts to collapse. Density and temperature 
of the very compressed material rise rapidly during the fast implosion of the core. 
When temperature is high enough the heavy nuclei synthesized during the whole 
stellar life start to be dissociated through photon processes: 

5'5Fe + 7 13^He + 4n- 124MeV, .^ 
^He + 7^2p + 2n-27.3MeV. ^ ' 

This endothermic photodisintegration processes start to drain energy from the core 
and trigger the collapse. Additionally electron capture on proton produces electron 
neutrinos and reduces the electron degeneracy pressure: 

p + e ^ n + Ue- (1.2) 

Collapse is not forever. When the inner part of the stellar core reaches nuclear 
densities p > 3 — 8 x 10^^ g/cm^ it suddenly stops. Not all the core feels the 



1 



2 



Explosion 




200 400 600 

Time after core bounce (msec) 



Figure 1.1: Schematic illustration of the collapse of a massive star core and the onset of 
the explosion, superposed to a diagram of the evolution of mass-shells from a supernova 
simulation 0/. The question-mark in the region between the neutrino-sphere (neutrino 
last scattering surface) and the stalled shock wave during the accretion phase indicates our 
present ignorance about the exact mechanisms that concur to produce the revival of the 
explosion. 



interruption (collapse velocities above the inner core are faster than acoustic waves) 
and matter continues infalling. In less than 1 msec a hydrodynamic shock wave 
forms and starts to move outwards sweeping the falling matter and reversing the 
collapse: the core bounces (see fig. Il.ip . 

At the very high densities of the collapsing stellar core neutrinos are trapped in the 
first hundreds kilometers where production and absorption processes are very fast. 
After a few milliseconds from the core bounce the shock wave reaches the region 
where neutrinos stream out freely. Electron neutrino luminosity suddenly increases 
by around a factor of 10 in less than 1 msec and the fluxes of z/^,r, ^/x,r are turned 
on (see sect. 12. ij) . Neutrino "fireworks" start. 

In tens of seconds more than 10^^ neutrinos of all flavors with energies in the range 
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10-50 MeV, are emitted carrying out almost all the gravitational binding energy 
released in the core collapse: 



where i?NS and Mns is the final mass of the neutron star formed in the process. 

In the inner regions of the massive star the future of the shock wave is uncertain. 
More than thirty years of core collapse simulations show that before leaving the 
core, the shock wave looses completely its energy through photodisintegration of 
nuclear matter, and becomes an almost static accretion front. Then what produces 
the supernova explosion? No complete consensus still exists in the answer of this 
important question [U although for sure, as was first realized twenty years ago 
j7] , neutrinos do play a central role. The absorption in the stalled accretion shock of 
even a small fraction of the huge neutrino fiux coming from the inner regions can be 
enough to revive the shock wave and restart the explosion, hundred of milliseconds 
after the initial core bounce. 

This delayed explosion scenario is by now the preferred way to explain to the "su- 
pernova explosion enigma" . However self consistent hydrodynamic simulations have 
not yet been able to obtain a successful explosion [HElini- The problem does not 
seem to be the mechanism in itself, in the sense that it is generally accepted that the 
delayed explosion occurs, and that neutrinos play a central role. More likely, other 
concomitant processes (convection, rotation, magnetic fields, etc.) could give the 
"last kick" . However the processes at the moment are not completely understood. 

In summary, neutrinos are copiously produced in supernova and (together with other 
not well understood concomitant processes) could finally trigger the explosion itself. 
However supernova neutrinos can actually do more than producing the explosion. 
The detection of these neutrinos can open a door to understand the physics of 
supernova by giving us an unique "snapshot" of the core collapse processes in a way 
that no other signal (except perhaps the still undetected gravitational radiation) 
can do. 

Supernova neutrinos are also very interesting for particle physics. On one hand 
they travel along huge distances, even in the astronomical sense (10 kpc is a typical 
number which is ~ 10^ times larger than the Earth-Sun distance), and any small 
neutrino mass will produce a tiny (but detectable) time delay among neutrinos of 
very different energies. On the other hand, its energy is in the right range to make 
neutrino oscillations across the Earth "visible" , allowing to measure neutrino mixing 
parameters with unprecedented sensitivity. 





CHAPTER 1. INTRODUCTION 



4 



In this work we will concentrate on the supernova neutrino potential to give us 
information on the absolute scale of neutrino masses. 

During the past few years, atmospheric [111311011111112] and solar [H } ITi l ITHl ITHl IT7| 

ITHl IT^ 1201 1^ neutrino experiments provided strong evidences for neutrino flavor 
oscillations and therefore for non vanishing neutrino masses. The KamLAND results 
\22\ I2S1 on the depletion of the Ve flux from nuclear power plants in Japan, and the 
K2K indication of a reduction in the flux from the KEK accelerator [2^, gave a 
final confirmation of this picture. 

However, to date all the evidences for neutrino masses come from oscillation ex- 
periments, that are only sensitive to squared mass differences and cannot give any 
information on the absolute scale of these masses. The challenge of measuring the 
absolute value of neutrino masses is presently being addressed by means of a remark- 
ably large number of different approaches, ranging from laboratory experiments to 
a plethora of methods that relay on astrophysical and cosmological considerations 
(for recent reviews see |25| |2HI ) . 

From the study of the end-point of the electron spectrum in tritium /5-decay, labora- 
tory experiments have set the limit m^^ < 2.2 eV j2Zll2Hl- If neutrinos are Majorana 
particles, the non observation of neutrino-less double /3 decay can constrain a partic- 
ular combination of the three neutrino masses. Interpretation of these experimental 
results is affected by theoretical uncertainties related to nuclear matrix elements cal- 
culations. This is reflected in some model dependence of the corresponding limits, 
that lie in the range mf < 0.2 - 1.3 eV 123 EOl 12111201 • 

Tight bounds Yli^i^i < 0-6 ~ 1.8 eV have been recently set using WMAP obser- 
vations of cosmic microwave background anisotropies, galaxies redshift surveys and 
other cosmological data (for a recent review see |221 and references therein). How- 
ever, these limits become much looser if the set of assumptions on which they rely is 
relaxed (see [H31IS11ES] for discussions on this point), and might even be completely 
evaded in exotic scenarios where neutrinos can annihilate into hypothetical light 
bosons, implying a suppression of their contribution to the cosmic matter density 
and negligible effects on structure formation at large scales |36j . 

Zatsepin j^T] was probably the first to realize that information on a neutrino mass 
could be provided by the detection of neutrinos from a Supernova explosion. Other 
early (independent) contributions in studying the potential of supernova were put 
forth in refs. [3H1 EH EO] • The basic idea relies on the time-of-fiight (TOF) delay 
Attof that a neutrino of mass m,^ and energy i?^ traveling a distance L would suffer 
with respect to a massless particle: 
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1 Tn,, 

Indeed, already in the past, the detection of about two dozens of neutrinos from 
SN1987 11111121 US] allowed to set upper limits on m^. Due to the low statistics, the 
model independent bounds derived were only at the level of m-p^ < 30 eV [Hj while 
more stringent hmits could be obtained only under specific assumptions jlHl IIEI 113 
HHj More recently, a detailed reexamination of the SN1987 neutrino signal based 
on a rigorous statistical analysis of the sparse data and on a Bayesian treatment 
of prior informations on the SN explosion mechanism, yielded the tighter bound 
mv, < 5.7 eV jlUj. 

The first observation of neutrinos from a SN triggered in the years following 1987 an 
intense research work aimed to refine the methods for neutrino mass measurements, 
in view of a future explosion within our Galaxy. With respect to SN1987, the 
time delay of neutrinos from a Galactic SN would be reduced by a factor of a few 
due to the shorter SN-earth distance. However, the neutrino flux on Earth would 
increase as the square of this factor and, most importantly, the large volumes of the 
neutrino detectors presently in operation will yield a huge gain in statistics. In recent 
years several proposal have been put forth to identify the best ways to measure the 
neutrino time-of-flight delays, given the present experimental facilities. Often, these 
approaches rely on the identification of "timing" events that are used as benchmarks 
for measuring the neutrino delays, as for example the emission of gravitational waves 
in coincidence with the neutrino burst [2010]; the short duration z/g neutronization 
peak that could allow to identify time smearing effects [31] , the sudden steep raise of 
the neutrino luminosity due to neutrino-sphere shock- wave breakout |H2] , the abrupt 
interruption of the neutrino flux due to a further collapse of the star core into a black 
hole [Sni El] )• The more robust and less model dependent limits achievable with 
these methods are at the level of mj, < 3 eV. Tighter limits are obtained only under 
specific assumptions for the original profiles of the SN neutrino emission, or for the 
astrophysical mechanisms that give rise to the benchmarks events. 

In this thesis we present a new approach to the measurement of the absolute scale of 
neutrino mass from a future Galactic supernova signal. The new method proposed, 
which is formally based on the principles of Bayesian reasoning, is rather independent 
on specific astrophysical assumptions and makes use of the whole statistics of the 
signal. The method has been tested under a wide range of different conditions for the 
neutrino emission, mixing schemes and for different present and planned neutrino 
detectors. The detailed results of such tests are presented in this dissertation and 
allow us to conclude that even with an approximate description of the neutrino flux 
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and spectrum, competitive mass limits can still be obtained from the analysis of the 
next Galactic supernova signal. 

We have organized this dissertation in four main chapters. In chapter |21 we review 
the information available about the expected characteristics of the emitted neutrino 
fluxes and spectra. We have put special attention on those characteristics that can 
be identified as robust predictions regarding the supernova neutrino emission. We 
discuss in some detail how the physical processes involved in the core collapse, shock 
wave emergence and revival determine these robust features. Chapter Elis dedicated 
to "construct" the detected signal according to the informations acquired in chap- 
ter 121 The effect of neutrino oscillations are also taken into account in this step. 
Neutrino flavor oscillations are determinant to model in a correct way an observed 
supernova signal. Although our aim will not be that of studying oscillations signa- 
tures in supernova, the description of this phenomenon is an obligatory aspect of 
any description of a neutrino signal. Using these informations we describe the pro- 
cedure to construct a neutrino signal, and we evaluate the expected main properties 
of the signal under several mixing and detection conditions. Chapter |3] is devoted 
to explain with more detail the different methods that have been devised to mea- 
sure or constrain a neutrino mass studying supernova neutrino. More importantly, 
we describe in this chapter with full detail the statistical method that we have put 
forth to measure a neutrino mass from the analysis of the whole statistics in a fu- 
ture supernova signal. The method has been tested with the procedures described 
in the final chapter |S[ where the results obtained from such tests are presented in 
full detail. A comprehensive comparison of various results and an analysis of the 
conclusions that can be drawn from our study is also contained in this last chapter. 

In Appendix ^ we describe the main ideas and mathematical formalism of the 
Bayesian approach to data analysis. Appendix ^ is devoted to present a simple 
computer tool (publicly available) that has been designed specifically to generate 
synthetic supernova signals. Some of the underlying algorithms and numerical tech- 
niques used in the development of our computer codes are also presented in this 
appendix. 

A complete list of symbols and abbreviations used along the text is given in Ap- 
pendix o 
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Neutrino fluxes and spectra 



Before proceeding to construct numerical simulations of a Supernova neutrino signal 
it is necessary to establish which are the features of the neutrino emitted signal for 
which the theory is reliable and which features are more uncertain. 

Although different Supernova simulations produce neutrino signals with different 
sets of properties, there are a set of overall features common to almost every result: 

1. The neutrino luminosities and fluxes evolve with time in a rather well un- 
derstood way. In the first tens of milliseconds the luminosity and flux of all 
neutrino flavors increase rapidly. After this, the luminosity and flux stabilize 
in their maxima and start to decay slowly for several hundred of milliseconds. 
A slow decay on a time-scale of several seconds follows and finally the signal 
turns off. 

2. The supernova core is the only place besides the early Universe where neutrinos 
are in thermal equilibrium. The observational consequence of this fact is the 
emission of a quasi-thermal neutrino spectrum. 

3. The average energies of different neutrino flavors are hierarchical, i.e. muon- 
and tau-(anti)neutrinos (hereafter we will refer them globally as u^) have a 
harder spectrum than the electron antineutrinos (z/g) which in turn are hotter 
than the electron neutrinos (z/g). 

Despite the wide range of detailed results obtained in supernova simulations during 
the last three decades, we can be confident that a future supernova signal will exhibit 
these properties. 
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In this chapter we will describe in some detail these robust features of the neutrino 
emission from supernovae, its physical motivation and how they are supported by 
self-consistent supernova simulations. In order to be prepared for the chapters where 
this information will be used to generate and analyze realistic supernova signals, 
we have attempted to translate some of the described properties into analytical 
expressions. 

A more comprehensive and complete description of the emitted neutrino properties 
can be found in refs. [521 EHl EZ] and references therein. 

2.1 Evolution of Neutrino fluxes 

One of the best known features of neutrino emission from Supernovae is the way 
the neutrino fluxes evolve with time (we are interested only in times after the core- 
bounce). 

Although the detailed time structure and some of the time-scales of the neutrino 
fluxes differ from one simulation to another, it is possible to recognize the existence 
of mainly three different phases in the evolution of the rates. Each of these phases 
are directly related to well known physical processes in the collapsed core of the 
star. 

The main phases of neutrino's fluxes evolution is illustrated in figure ITTl 
2.1.1 Fast rising - the shock breakout 

Once the inner region of the collapsing stellar core reach supra-nuclear densities a 
hydrodynamical shock wave forms and starts to travel out. After less than few mil- 
liseconds the shock wave reach the region of the core where neutrinos escape freely. 
A sudden increase in the plasma temperature produced by the energy deposited by 
the shock, together with the sudden transition from coherent scattering on "iron" 
nuclei to scattering on nucleons produce a huge release of electron neutrinos (the 
prompt neutrino burst oi deleptonization burst). 

The shock breakout turn on the emission of in less than 1 msec (neutrino-sphere 
crossing time of the shock wave [Sn]). After this turn-on the heat from the shocked 
material above the neutrino-sphere and below the shock wave produces a steadily 
increase in the flux that last from 20 to 50 msec (see figure 12. 1|) 
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Figure 2.1: Neutrino fluxes evolution could he described with three main phases: 
shock breakout (fast rising), accretion and mantle cooling (maximum and plateau) and 
Kelvin- Helmholtz cooling (decay). The fluxes were taken from simulations by Totani et al. 



2.1.2 Flux maximum - accretion, mantle cooling and explosion 

Neutrino emission and nuclear dissociation consume energy of the emergent shock 
wave. Even before it leaves the iron core, the shock wave becomes an almost sta- 
tionary accretion front. In these conditions a new source of neutrinos appears. The 
inner core (protoneutron star) emits its heat through neutrinos in a steady fashion. 
The same thing happens to the hot material below the accretion front (core mantle). 
However, as new material is accreted in the stalled shock, it is dissociated and pro- 
duce extra by electron capture on protons. The hot material in the shock front 
also produces electron-positron pairs and the capture of positrons on neutrons, in 
the neutron rich material of the mantle, creates an additional flux. The emission 
of neutrinos during this phase is characterized by a slow decay in the flux on a 
time-scale of one hundred of milliseconds. 

The emission of neutrinos by accretion is suddenly truncated when the shock wave 
revives and the delayed explosion starts. 
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2.1.3 Decay - Kelvin-Helmholtz cooling 

Once the revitalized shock wave emerges a more stable and long phase starts. The 
gravitational binding energy of the newly formed protoneutron star ()1.3|) is radi- 
ated as neutrinos on time-scales comparable with the diffusion time-scale {Kelvin- 
Helmholtz cooling). 

This phase is characterized by a power law decrease in the neutrino flux and the 
emission of an important fraction of the protoneutron star binding energy. The 
Kelvin-Helmholtz cooling phase last for tens of seconds to minutes and it ends when 
the material of the newly formed neutron star becomes transparent to neutrinos. 

As was mentioned before this basic description of the overall features of the emitted 
neutrino flux is a robust prediction of supernova models. Its global features have 
been confirmed by the observation of the neutrino signal from SN1987A. The initial 
short burst of neutrinos (~ 1 sec) observed in the signal is consistent with an ac- 
cretion phase and the observed total duration of the signal (~ 10 sec) confirms 
the emission of neutrinos from a cooling neutron star [HH] . 

2.1.4 Analytical flux model 

In order to quantify the overall features of the flux described above it is possible 
to construct a simplified model for the neutrino emission. Such a model will be an 
useful tool to understand the characteristics of fluxes obtained in simulations. In 
section I4.5.H| we will start with this analytical approximation to construct a test 
model to analyze a realistic neutrino signal. 

We can recognize two main components in neutrino emission: 

• Cooling component. On a crude level of approximation neutrino emission 
from the supernova core and mantle could be written in terms of the Stefan- 
Bolt zmann law (see (SHI for a detailed discussion), 



where -Refr is the effective radius where neutrinos are emitted and a quasi- 
thermal spectral distribution F{E] Tcfj (t)) is assumed, with Teg{t) as the spec- 
tral temperature. Cc is a constant in which we have absorbed geometrical and 
physical terms. 



dn 
Itt 



em 



47ri?eff(t)' / F{E-%s{t))dE = C,R,s{tf%s{t)\ 



(2.1) 
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• Accretion component. When accretion starts an additional contribution to the 
flux is turned on. /5-processes producing and neutrinos in the accreting 
matter are in local thermal equilibrium (LTE) and thus neutrino spectra could 
be also described as thermal with a temperature equal to the accreted material 
temperature T^. The neutrino emission rate per unit of mass of accreted 
material will be given by |49j . 



The E"^ term arises from the energy dependence of the /^-processes cross- 
section. Ca{Ta) is a parameter depending on the spectral temperature of the 
accreted matter that can be assumed constant in time. Ma{t) is the mass of 
accreted material at time t. 

When the cooling emission component 1^2.1^ and the accretion component ()2.2j) are 
summed up the total flux on Ue and Ve is 



As can be seen the temporal evolution of neutrino flux is a function of the temporal 
evolution of the effective emission radius Res, the effective thermal temperature T^s 
and the rate of matter accretion M^. Different phenomenological models could be 
formulated to describe the behavior of this quantities. For the cooling component 
it is common to use a exponential cooling law, i.e. Teg = T^e"*/^^ and assume an 
almost constant effective neutrino-sphere radius. This kind of models seems to flt 
well the gross behavior of the SN1987A signal [SH]. In a more complete and detailed 
reanalysis of SN1987A signal Loredo and Lamb [IH] has recently found that a power- 
law cooling model, T{t) = To(l -|- t/4'~fT)~'^ which is capable to reproduce a wide 
range of cooling behaviors (e.g. when 7 — >■ oo we obtain an exponential cooling) 
could flt even better the signal jlH]. 

For the accretion component a truncated power law decay for the mass of accreted 
matter, M„(t) = M^e-^*/*")"'' (1 + t/tfe)" seems to flt well supernova simulations 
results HH]. 

Using the preferred phenomenological models the flnal neutrino flux for times pos- 
terior to the shock breakout ()2.3p is written as 




(2.2) 



dn' 



C,R,s{tf%s{tf + Ca{Ta)Ma{t). 



(2.3) 



dt 



dn' 



C 



e(-*/ta)"- 



(2.4) 



dt 



(l + t/7r)T (l + t/4) 
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where C = CcTo and A = Ca{Ta)Mo and we have reparametrized the cooling com- 
ponent in terms of ric = 47 and tc = 7t. Note that this flux does not behaves 
properly for t ^ since it is aimed to describe the emission of neutrinos for times 
after the shock breakout. In order to use this analytical approximation to flt a signal 
we will modify ()2.3p in section 14.5.31 to ensure the proper asymptotic behavior of 
this function. 

In figure 17!^ we compare the and fluxes from supernova simulations with this 
analytical model. As can be seen the flux model ()2.4|) reproduces well the behavior 
of the flux in the accretion and cooling phases. 



3e+57 



2.5e+57 



Buras et al. 2004 
Woosley et al. 1994 
Analytical fit 




2e+57 



1.5e+57 



le+57 - 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

t (sec) 

Figure 2.2: Schematic analytical fits ofVe fluxes from two supernova simulations using 
the flux model \2.4\l - 



Finally in flgure 12.31 we present fluxes of all neutrino flavors coming from three 
different supernova simulations. 
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5e + 56 i' 




Figure 2.3: Supernova neutrino fluxes from three different self- consistent simulations: 
Totani et al. \5<^ . Woosley et al. 16 (J^ and Buras et al. l&lf - The main overall features of 
the fluxes described in the text could be recognized here. 
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2.2 Neutrino Spectra 

Neutrinos trapped inside the inner regions of a supernova core are maintained in 
thermal equihbrium with the surrounding plasma through inelastic scattering. For 
this reason it is expected that they will be emitted with a quasi-thermal spectrum 
characterized by a spectral temperature close to the temperature of the plasma in 
the region where neutrinos decouple energetically. More than 30 years of analytical 
and numerical studies of neutrino spectra formation in supernova have confirmed 
this basic picture. However, the fine details of neutrino transport and emission near 
to the last scattering surface are highly non trivial and the precise determination of 
the spectral properties of the emitted neutrinos is not a simple task. 

Let us examine here the main features expected for the neutrino spectra and its 
origin. A more comprehensive and complete discussion could be found in ref. jHZj 
and references therein. 

2.2.1 Electron neutrinos and antineutrinos 

The production and transport of and Ue is simple when compared with the trans- 
port of Vx- Due to the presence of electron and positrons, Ve and Ve can undergo 
charge current reactions (/5-processes ep ^ nue, e+ra ^ pUe) that are almost ab- 
sent for all other flavors. With a cross-section much larger than other scattering 
processes, /3-processes dominate completely the transport of Ue and I/g inside the su- 
pernova core. The continuous creation and capture of neutrinos in these processes 
inside the core maintain and Ve in LTE with the surrounding medium. 

Neutrinos diffuse outward until the /5-processes become ineffective at the neutrino- 
sphere. Since the cross section of /3-processes are energy dependent the neutrino- 
sphere radius and therefore the characteristic thermal temperature for i>e and Vg, 
are different for neutrinos of different energies. High energy neutrinos will stream 
out from larger radius with lower temperatures while low energy neutrinos will 
be radiated from a lower radius and higher temperatures. Therefore the energy 
distribution of the emitted neutrinos will not be exactly thermal. The emitted 
neutrino spectra will be "pinched", i.e. high energy and low energy tails of the 
spectrum will be depleted due to the different thermal temperatures characteristic 
of the radius where neutrinos of different energies are emitted (see fig. 12. 4^ . 
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Spectrum at 100 msec ^ 




10 20 30 40 50 60 

E(MeV) 

Figure 2.4: I7g spectrum at 100 msec (accretion and mantle cooling phase) taken from the 
simulation of Buras et al. Wlf (histogram). The dotted curve correspond to a Fermi-Dirac 
fit with zero chemical potential. The pinching effect, i.e. the depletion of the distribu- 
tion at high and low energy, can he clearly recognized in the distorted Fermi-Dirac and 
a- distribution fits. As can be seen, for this particular case the a- distribution fits slightly 
better the spectrum than the pinched Fermi-Dirac. 



2.2.2 Non-electron neutrinos 

The emission of from tlie supernova requires a more complicated description. In 
this case different neutral current reactions play comparable roles in the neutrino 
transport. 

Vx can undergo three kind of processes: creation reactions, energy exchange scatter- 
ing and isoenergetic scattering. 

Deep inside the core neutrinos are created and annihilated through three main 
processes: nuclear bremsstrahlung NN NNuu, electron-positron annihilation 
ee~^ ^ VxVx and neutrino-neutrino annihilation Vei^e ^ ^x'^x- There, neutrinos are 
in LTE and trapped. At a given radius creation/annihilation processes freeze out 
and no more neutrinos are created. This condition defines what is called the number 
sphere. The emergent neutrino total flux is fixed in this layer. Outside the number 
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sphere neutrinos still undergo energy exchange scattering on electrons, positrons 
and nucleons and maintain some thermal contact with the plasma. When density 
and temperature is lower inelastic scattering on nucleons and becomes inefficient 
and neutrinos decouple energetically from the surrounding medium. This condition 
defines the so-called energy sphere. Above the energy sphere neutrinos still interact 
strongly with the plasma through elastic scattering on nucleons. Finally at some 
radius elastic scattering becomes in turn inefficient and neutrinos stream out from 
the neutrino-sphere or transport sphere. In figure we illustrate schematically this 
process. 

What will be the emergent energy distribution of z/j,? Naively one could think that 
since thermal contact is maintained only inside the energy-sphere neutrino average 
energy will be related with the plasma temperature at that layer. However the 
energy dependence of the nucleon scattering outside the number-sphere acts as a 
filter jn^l lowering the average energy to 50-60% its value close to the energy sphere. 

The complex interplay of all factors which determine the emergent fiux makes 
harder to predict exactly which will be the spectral properties of these fiavors. We 
will come back on this issue in section IT^ 



2.2.3 Spectral pinching 

The quasi-thermal nature of the neutrino spectra is a robust feature of neutrino 
emission from a supernova. 

There are two ways to quantify the non-thermal distortion (pinching effect) of the 
neutrino spectra. On one hand it is possible to use higher order momenta of the 
neutrino energy distribution that contain informations about the high energy tail. 
Raffelt [02] first quantifies the pinching effect by introducing a "pinching parameter" 
defined in terms of the ratio of the second and first momentum 



a (E) 

Where a ~ 1.3029 is a weight constant equal to the ratio of the second and first 
momentum of a Fermi-Dirac distribution with zero chemical potential. In the ab- 
sence of any pinching and without degeneracy p = 1. A pinched spectrum will be 
characterized by p < 1 while an anti-pinching, i.e. an enhancement of the spectrum 
tails, corresponds to p > 1. 

On the other hand has been common in the literature to fit the supernova neutrino 
spectrum using a Fermi-Dirac distribution with an effective degeneracy parameter 
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Figure 2.5: Schematic representation of transport. Neutrinos are produced and an- 
nihilate through nucleon hremssthralung and pair processes. When creation/ annihilation 
freeze-out at the number sphere the total flux is fixed (number of arrows). Neutrinos 
maintain a thermal contact with the medium until they reach the energy sphere (style of 
lines indicates the neutrino average energy). Finally above the transport sphere (neu- 
trino-sphere) neutrinos freely stream out. Its average energy at this layer (style of arrows) 
is close to that at the energy sphere. 
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rj that enters the distribution like an effective chemical potential: 



FFD(i?;T,ry)oc^^^^. (2.6) 

The introduction of rj allows to fit the distorted tails of the distribution. Its value 
could therefore used as a measure of the spectrum pinching. 

Recently an alternative two-parameter distribution has been proposed ES] '■ 

[E; e, a) = N{e, a) {E/ef 6-^"+^) , (2.7) 

where the normalization constant is given by A^(e, a) = (a + l)"+^/r(a + l)e. It 
has been shown that this new distribution, denoted a- distribution, in j^j fits slightly 
better the neutrino spectra and has the nice property of allowing a simple analytical 
estimation of the two spectral parameters e and a directly in terms of the first 
and second momentum of the energy distribution. Using the well known relation 
aV^a) = T{a + 1) it is easy to verify that: 



2 + a _ {E^ 
l + a~ {Ey 

Using these relations, the pinching parameter p could be expressed in terms of a as: 



ap= (2 + «)/(! + a). (2.9) 

Therefore when an a-distribution is used to fit the neutrino spectrum the value of 
the a parameter can quantify the spectrum pinching. 

Figure 12.41 illustrate the quality of fits of a numerical neutrino spectrum with a 
Fermi-Dirac and an a-distribution. 

The relation between p, rj and a is depicted in figure 12.61 Typical values for a from 
a self-consistent supernova simulation are presented in figure 12.71 As can be seen 
the spectral pinching evolves to lower values as the protoneutron star and the core 
mantle settles and neutrino emission is confined to more compact regions jS]. 

Typical ranges for the spectral pinching obtained in simulations are [BIJ 0] : 



p,^ ^ 0.88 0.92 , r/,^ ^ 3 ^ 5 , a,, ~ 4.0 ^ 5.6 

~ 0.90 - 0.95 , 7/^^-2-4 , ^3.0-5.0 (2.10) 
p,^ -0.95 -1.00 , r/,^~0^2 , a,^ ~ 2.3 - 3.2 
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Figure 2.7: Evolution of the a parameter in the supernova simulation by Buras et al. I6'j) /. 
The pinching gets reduced at later times when neutrinos are emitted from more compact 
regions inside the protoneutron star. 
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2.3 Hierarchy of average energies 

The last robust feature of the emission of neutrinos from supernova is the hierarchy 
of the average energies of different neutrino flavors. 

As was discussed in the previous section the spectral properties of neutrinos emitted 
from the supernova core are determined by the transport processes near the neutrino- 
spheres. In the simple case of v^. and z7e the average energy is determined by the 
temperature of the plasma around the neutrino-spheres. Although the cross section 
for /3-processes (the main source of opacity for these flavors) has the same energy 
dependence for v^. and T/g, the asymmetry between protons and neutrons in the 
supernova core produces a significant difference in the opacity of antineutrinos and 
neutrinos. Since neutrons are more abundant, z/g decouple from matter at a larger 
radius than z/g, and therefore has a lower average energy, i.e. {E^^ < {E-pJ. 

The size of the difference between the average energies of z/g and z/g depends basi- 
cally on the core structure (density, temperature and lepton number profile) which 
determines the position of the respective neutrino-spheres. In supernova simulation 
typical values ranges from 30% jHHEn] to 50% jHOJE^] (see figure . 

On the other hand the transport of is dominated by neutral current reactions that 
have smaller rates than the /9-processes. This implies that decouples at a lower 
radius than Ve and therefore it is expected that {E^J > {E-pJ > {E^J. However 
the determination of how large are the differences between the average energies of 
and Vf, is not easy. If the average energies were only determined by the spectral 
temperature at the energy-sphere, large differences would be expected. However 
the filter effect of the elastic reactions outside the energy-sphere for contribute 
reducing the average energy. It has been shown in ref. jH] that the final difference 
will depend on the details of the transport in that regions. 

Most of the self-consistent simulations have made use of simplified descriptions of 
i^x transport mainly due to the fact that does not play an important role in 
the supernova dynamics. In these simulations the resulting spectra has typical 
energies differences that range from 30% up to 100%! (see figures l?!Hk .b). 

Large spectral differences of this size have been exploited in several works to study 
the potential information contained on a supernova signal. Minakata et al. [Hi] 
have claimed that an inverse hierarchy of neutrino masses could be rejected using 
the informations contained in the SN1987A signal. Lunardini and Smirnov jHEl have 
derived similar conclusions about the potential of SN1987A signal to give informa- 
tion about neutrino oscillations parameters. The potential of a future supernova 
signal has also be studied in the presence of large spectral differences |67j . 
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Figure 2.8: Evolution of the mean energy from three different supernova simulations. It is apparent the big difference 
between the "standard" result of large spectral differences in simulations where a simplified description of transport 
was used (panels a and b) compared with almost equal average energies of Vx andV^, when a more precise and complete 
description of interactions are included (panel c). 
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Recently the inclusion of a detailed description of Ux microphysics in self- consistent 
hydrodynamical supernova simulations (see ref. [HZI has suggested that the stan- 
dard picture of large spectral differences could be wrong. Including all the relevant 
reactions in the transport, and other effects that had not been taken into ac- 
count in previous simulations, the average energy results almost equal to that of 
Ve- Although a partial hierarchy could still be present (at least during part of the 
emission) differences will not be larger than 20%, with 10% as a typical value j57j . 
Panel c in figure ITHl illustrates this point. 

We conclude that, according to the most recent results, any attempt to obtain in- 
formations on neutrino mixings and their mass hierarchy from a Galactic supernova 
signal will have to take into account that the differences between the neutrino ener- 
gies will be probably not very large. 
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supernova 



Using what we have learned about neutrino emission in Supernova we can now 
proceed to calculate the expected characteristics of a future supernova signal. One 
main part of this work regards using this description to create artificial signals that 
will be later used to test the statistical method that will be presented in the next 
chapter. 

In this chapter we present the theoretical description of a supernova neutrino signal 
and the general strategies that can be used to generate synthetic realizations of such 
signals. In section ITT] we present the most relevant ideas about supernova neutrino 
oscillations which will play a central role in the determination of the final features 
of any neutrino signal. Section 13.21 introduces the formalism and describes in some 
details how we model the neutrino detection process. In section 13.31 we compute 
the most relevant features of a future supernova neutrino signal, assuming different 
emission models and oscillation effects and analyze the interplay of the different 
aspects which are involved in the determination of the signal characteristics. Finally 
section 13.41 describes the particular strategies, algorithms and numerical techniques 
used to create synthetic signals. 

The material presented in this chapter is a collection of useful results obtained from 
rather standard descriptions of neutrino propagation and detection phenomenology. 
We intend just to present synthetically the theoretical results behind the computer 
codes used in our work. 
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3.1 Oscillations of supernova neutrinos 

Neutrinos will undergo oscillations on their way from the supernova core to a de- 
tector. First, when they traverse the outer shells of the stellar core and then when 
they arrive to the Earth and travel through our planet's mantle and/or core before 
being detected. 

In the last two decades many works has been devoted to study supernova neutrino 
oscillation phenomenology. Many of these efforts were devoted to understand how 
neutrino oscillations affected the SN1987A signal jHSHHn]- More recently the effect 
of neutrino oscillations on a supernova signal have been studied in the attempt 
to determine how the distortions caused by oscillations could be used to estimate 
and / or constrain the neutrino mixing parameters [HZl EHl EHl IZOl 1^ • 

Here we will summarize the most relevant facts and predictions of these works. More 
detailed discussions can be found fro example in refs. [HZIinn] and references therein. 

3.1.1 Neutrino Oscillations - Basic facts 

Neutrinos are produced and detected as weak-interaction (flavor) eigenstates Va {ot = 
e,fj,,T). We know that the flavor eigenstates do not diagonalize the mass matrix, 
and after being produced they propagate as a superposition of mass eigenstates 
Ui {i = 1,2,3). The flavor and mass neutrino basis are related through a unitary 
transformation: 



where we denote i'w= i^^e^n^^T)"^ , ^m= {.^i^2^z)^ and U is the Maki-Nakagawa- 
Sakata matrix (MNS). U (whose elements we will denote as Uai) can be param- 
eterized in terms of three mixing angles ^12, ^23 and ^13, and assuming 5 = it 
reads [Za US]: 



where we have used the notation Cij = cos% and Sij = sin 6ij. 

The propagation of an ultra relativistic neutrino mass eigenstate E ^ p{l +m1/2E) 
obeys an approximate Schrodinger-like equation of mot ion [7^: 



(3.1) 




(3.2) 
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I— ~ Hmi^m, (3.3) 

where the Hamiltonian in the mass basis is Hm = (l/2£')diag (m^, m^, mg). In 
the flavor basis the propagation equation could be obtain from ()3.3|) applying the 
transformation ()3.1|) . 



= HMiU^^y^, (3.4) 



dt 

. duw 
I —— 
dt 



Hwi^w (3.5) 



with the Hamiltonian defined as Hw = U Hm ■ 

The evolution of a mass eigenstate is given by z/j(t) = exp[—{m]/2E)t]i'i{0) and 
can be used to write down how a neutrino which is produced as a definite fiavor 
eigenstate i^o(O) = J^^aiPiiO) evolves: 

'ya{t)=J2Uaie~'^'iy,{0). (3.6) 

We see that the neutrino fiavor (given by the particular mix of mass eigenstates) 
changes along the path, provided the neutrino masses are different. This is the 
reason underlying neutrino vacuum oscillations. 

The projection of the neutrino state at t on the fiavor basis allows us to find the 
general conversion probability in vacuum: 



^^f, = P(u^ ^ Up) = \{upMt))f (3.7) 



\j2U(^jU*Mu,{0))\ (3.8) 
Y,U;,U^,Uf,,U:,e-^^\ (3.9) 



where Amfj = mf — m| are the squared mass differences. The above equation tells 
us that a neutrino fiavor oscillate with a spatial period £osc = 4:7rE/ Am'^j called the 
oscillation length. 

In matter, the effect of coherent neutrino interaction with the medium can be ac- 
counted for by adding a potential term to the equation of motion (|3.5|) . 
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w 



dt 





fy 









Hw+ 1 











1 









0) 





(3.10) 
(3.11) 



where V = \f2Gpne is the effective electron neutrino potential [—V for antineu- 
trinos) that is due to coherent charge current scattering off electrons. In the last 
equation we have redefined the fields to absorb the neutral current effective potential 
which is common to all flavors and contributes an overall phase. 

With the introduction of the interacting terms the vacuum mixing matrix U does not 
diagonalize the weak Hamiltonian and therefore the vacuum mass eigenstates are 
no longer propagation eigenstates. A new matrix t/™ which diagonalize must 
be introduced, by defining a new instantaneous mass eigenstates basis z^= W^i/^ 
which in general is different from point to point along the neutrino path if the density 
and therefore V changes. 

In the two flavor case (which will be interesting for supernova mantle oscillations 
described in section I3.1.3|) the effective weak Hamiltonian is diagonalized by 
the in-matter mixing matrix 



U 



cos 6m sin 6„ 
— sin 9,m cos 9r, 



(3.12) 



where 6m is the mixing angle in matter and reads: 

. sin 2^ 

sm 2u„ 



Vsin^ 26 + {cos2e - 2EV/Am^y ' 
The eigenvalues of U^H^U^'^ give the squared effective neutrino masses ^\ 



(3.13) 
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ml + ml 2EV Am^ 



mf + m2 



2 2 
2EV Am2 



sin^26 + {cos26 -2EV/Am^y , (3.14) 
sm^26 + {cos26 -2EV/Am'^)^. (3.15) 



From this we get Afi^ = Am^^/sm'^26 + (cos 2^ - 2EV/Am^y and O = AnE/Afx. 
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In summary, the effect of the different neutrino interaction with the surrounding 
medium is that of producing a change in the mixing angle and therefore in the 
oscillation rate This change is maximal when the resonant condition is fulfilled: 



2EV = Am^cos2e. (3.16) 

When the medium density changes another "degree of freedom" is added to the 
oscillation dynamics of neutrinos. Starting from the equation of motion for flavor 
eigenstates in matter ()3.11|) we can obtain the instantaneous mass eigenstates i^by 
solving: 



KU^^l^, (3.17) 
= 777;diag(/ii,/i2,/X3)-^t/'"^h^ (3.18) 



dt \2E 

The variation of along the neutrino path makes the Hamiltonian in the last 
equation not diagonal and therefore z^are not longer propagation eigenstates. As 
a result not only neutrino flavor change due to differences in the phase evolution of 
the propagation eigenstates as in vacuum, also the states mix among themselves. 

The mixing of instantaneous mass eigenstates is governed by the ratio between the 
diagonal and off-diagonal terms of the Hamiltonian in eq. ()3.18p . The off-diagonal 
terms which are provided by dV^^ j dt depend on the rate of variation of the in-matter 
mixing angle d9m/dt which is determined by the density gradient. 

Near the resonance region ()3.16p the diagonal terms of the effective mass Hamilto- 
nian are almost equal and cancel out by a phase redefinition of the fields. Therefore 
it is in that region where the instantaneous mass eigenstate are more important. 
If density varies slowly inside the resonance region (and the off-diagonal terms are 
close to zero) the instantaneous mass eigenstates are almost propagation eigenstates 
and do not oscillate. It is said that neutrinos perform an adiabatic transition. But, 
if the density gradient is steep enough, the instantaneous mass eigenstates strongly 
oscillates and neutrinos could emerge from the resonance layer as a different instan- 
taneous mass eigenstate. It is said that the neutrino state jumps. 

Quantitatively, the adiabaticity condition can be expressed in terms of the mixing 
parameters and the medium properties as |74j : 



_ Am^ sin^ 26 



d log Ue 



dt 



(3.19) 
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where 7 is called the adiabaticity parameter. 

The jumping probability of instantaneous mass eigenstates near the resonance region 
is approximately given by the Landau- Zener Formula^: 

Pj = exp(-7r7/2). (3.20) 

The resonant conversion of neutrinos in matter through the above process is con- 
ventionally called the Mikheev-Smirnov-Wolfenstein effect (MSW) |76[ I77j. 

3.1.2 Neutrino mixing parameters - experimental values and limits 

As we saw in the previous section the dynamics of active neutrino oscillation depends 
on 3 pairs of parameters: {Aml2,0i2) , (^'^23'^23), (^^135^13)- 

The study of neutrino oscillations in solar, atmospheric, reactor and accelerator 
experiments has allowed us to measure accurately some of these parameters, and to 
put limits on other ones. In the following, we present the most recent values and 
limits obtained from a global analysis of all the available experimental results |78j : 

• [Am?2 = 7.4 - 9.2 x 10~^eV\ Ou = 28° - 37°] (99.7% CL). 
This is the so-called LMA solution to the solar neutrino problem. 
It is common to identify this pair of parameters as (Am^, ^0) 

• [lAm^gl = 1.5 - 3.9 X 10-^eV\ 623 = 34° - 45°] (99.7% CL). 

The sign of Am"^^ is still unknown and will depend on what is the mass hi- 
erarchy: normal hierarchy (NH) Am^g > or inverted hierarchy (IH) 
Am23 < 0. 

It is common to identify this pair of parameters as (Am^^^^, 6'atm)- 

• [Amfg ~ Am^g, < 12.9°] (99.7% CL). 

3.1.3 Neutrino Oscillations in the Supernova Mantle 

In the following sections we will apply these basic results to study neutrino oscilla- 
tions in the supernova mantle. 

^In principle the Landau-Zener formula is accurate under the assumption of a small value of 
the adiabaticity parameter and a linear variation of the density inside the resonance region. We 
have verified that for typical supernova density profiles and mixing parameters this formula works 
as well as the double exponential formula used in other works jZH] • 
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3.1.4 Initial conditions 



Inside the supernova core where neutrinos are produced, the effective neutrino po- 
tential felt by the electron (anti)neutrinos overwhelms the off-diagonal terms of the 
effective weak Hamiltonian ()3.1H1 and fiavor eigenstates are approximately lo- 
cal mass eigenstate. The correspondence between flavor and mass eigenstates in the 
production region will depend on the mass hierarchy jHZj: 



• with normal hierarchy we will have: 



for inverted hierarchy: 



V 



(3.22) 



Where u'^, u'^, V'^ and V'^ are redefined states which diagonalize the (z/^, i/^) sub- 
matrix of jnZj. It must be noticed that such redefinition does not change the 
supernova physics since for all practical purposes and u^. are phenomenological 
identical in supernova conditions. 



3.1.4.1 Resonance regions 

Neutrinos emerging from the supernova core traverse matter with densities ranging 
from 10^ g/cw? to almost 0. In the process they will find the required conditions for 
resonant conversion. Two possible resonances are allowed by the mixing parameters: 

• a high density H-Resonance corresponding to enhanced conversion into the 
heaviest mass eigenstate z/3 and determined by the condition: 

2E V{p'^') = Aml^ cos 2^13. (3.23) 

• and a low density L-Resonance associated to transitions between the lightest 
eigenstates: 

2E V{pT) = Am^2 cos 2^12. (3.24) 
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Since Aml2 > the L-Resonance will be only in the neutrino channel. 

On the other hand the H- Resonance will affect neutrinos and antineutrinos according 
to the mass hierarchy: 

• for normal hierarchy (Am23 > 0) H-Resonance will just affect neutrinos, 

• while in the case of an inverted hierarchy (Amgs < 0) H-Resonance will 
occur in the antineutrino channel. 

3.1.4.2 Factorization of Dynamics 

The hierarchical nature of the squared- mass differences Am\.^ : Amfg — 100 : 1 will 
produce a natural decoupling of the conversion dynamics in the two resonance layers 
inside the supernova mantle. 

At high densities (H-Resonance) the in-matter mixing of the lightest eigenstates 
will be suppressed. Transitions will affect just the heaviest eigenstate and the 

conversion description is reduced to the two flavor case (2/1,2/3). At low densities 

(L-Resonance) the mixing of the third eigenstate is near to its value in vacuum 
^ Ue3 ^ 10"^ and almost constant. Only the first and second eigenstate will 

mix and again the problem is reduced to a two flavor case (1^1, 1^2)- 

Additionally as Ues is very small the total survival probabilities for z/g and Vf. can 
be factorized [731^: 

P^PhX Pl. (3.25) 

Factorization of dynamics is a property which simplifies considerably the computa- 
tion of the total conversion probabilities in the supernova mantle without a detailed 
numerical study of the neutrino propagation equation. 

3.1.4.3 Jumping probabilities 

We have computed generic jump probabilities using the Landau-Zener formula ()3.2Up 
for typical conditions found in the supernova mantle and generic mixing parame- 
ters. Contours of equal probability in the (Am^,sin^^) plane for different neutrino 
energies, are depicted in figure ITTl 

Observing the contours, we see that the neutrino transitions in the L-Resonance, 
which are governed by the solar mixing parameters (Am^j, ^12) (LMA), is completely 
adiabatic, i.e. Pl = 0. 
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The case for the H-Resonance is not as simple. Since the neutrino transition through 
that layer is governed by the unknown mixing angle 6*13, the adiabaticity will be also 
uncertain. 

We can recognize two extreme cases |67j : 

• sin^ 6*13 ^ 10^'' (6*13 ^ 1°). In this case for all the interesting neutrino energies 
the adiabaticity condition is fulfilled, i.e. Ph — 0. 

• sin^ 6*13 < 10~^ (6'i3 < 0.06°). In this case adiabaticity is strongly violated 
implying Ph ~ 1 for almost the entire range of supernova neutrino energies. 

For 10~^ ^ sin^ ^13 ^ 10~^ we will have an intermediate behavior with jump prob- 
abilities ranging between 0.1 and 0.9. 

The energy dependency of the jump probability in the resonance region is another 
matter of concern. We have plotted in figure IT^ Pw as a function of neutrino energy 
for different values of the mixing parameter sin^ 6*13. It can be seen that for the range 
of interesting neutrino energies the jump probability in the extreme cases is almost 
constant. In the intermediate range of values of the mixing parameter, Ph will have 
only a mild variation in the energy range where most of the neutrino events will be 
expected. 

The above conclusions apply also to antineutrinos in the inverted mass hierarchy 
case where the jumping probability Ph becomes identical to Ph- 

3.1.4.4 Conversion probabilities 

Using the jumping probabilities and the fact that the dynamics in the resonance 
region decouple we can compute the transition probability {pa.i)Paii i-e. the prob- 
ability that a given (anti) neutrino flavor a emerges from the supernova in the i^^ 
mass eigenstate. Figure illustrates schematically the way these probabilities are 
computed following the evolution of the neutrino state from the emission region to 
the star surface. 

In table EH) we summarize the results of applying rules as those presented in fig. 13.31 
to compute transition probabilities for neutrinos and antineutrinos, assuming normal 
and inverted hierarchies. 

Transition probabilities allow for a straightforward computation of the final conver- 
sion probabilities, i.e. the probability that a given neutrino flavor a is detected as 
flavor (3. In terms of Pai the conversion probabilities are simply given by 
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Table 3.1: Transition probabilities 
Normal mass hierarchy 



Neutrinos Antineutrinos 





Z/i U2 


1^3 






^2 






PhPl Ph{1-Pl) 


(1 - Ph) 




1 










(1 - Pl) Pl 










1 







{1-Ph)Pl {1-Ph){1-Pl) 


Ph 










1 



Inverted mass hierarchy 



Neutrinos Antineutrinos 













1^1 


1^2 


1^3 







Pl 


(1 - Pl) 


Ve 


Ph 





(1 - Ph) 







(1 - Pl) 


Pl 







1 







1 








Vr 


i^-Pn) 





Ph 



Paf3 = y^^VaiPip- (3.26) 

i 

Where pij3 is the probability that the mass eigenstate i^^ be detected as the neutrino 
of flavor (3. 

In the absence of any other mixing effect neutrinos emitted from the supernova 
arrive to the Earth surface as incoherent mass eigenstates and is the amplitude 
of the projection on the flavor basis in vacuum: 

P^P = (3.27) 

We will focus on electron antineutrinos for reasons that we will give later. 

As an example, the electron antineutrino conversion probabilities for normal and 
inverted hierarchy read: 



P— 

^ ee 



Pl 
P^ 



fie 



\Uel\ 



e2\ 
\Ue3\ _ 
PH\Uelf + {l-PH)\Ue3f 

\Ue2f _ 
(l-PH)\Uelf + PH\Ue3\' 



NH 



IH 



(3.28) 
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We can see that in the case of a normal hierarchy the observed I^g will be an 
admixture of neutrinos emitted originally as V^. (Pge = cos^9i2 — 0.26) and V'^ 
{Pjce = cos'^612 — 0.74). On the other hand in the inverted hierarchy case the 
mixing will mainly depend on the adiabaticity of antineutrino transitions in the H- 
Resonance. If transitions are adiabatic, Ph — 0, the survival probability for Ve will 
be suppressed Pee — If^esl^ ~ 10~^ and the observed Ve events will be and Vr 
neutrinos converted in the supernova mantle into Ve {Pjm + Pre — 1). 

Similar considerations can be done for electron neutrinos (see ref. |67j). 
3.1.5 Neutrino oscillations in the Earth 

If neutrinos travel inside the Earth oscillations will mix the arriving mass eigenstates 
and the probabilities pij3 will be more complex. 

Due to the hierarchical nature of the mass spectrum, antineutrino oscillations in the 
Earth could be also reduced to a 2 flavor oscillation problem [HZI- In this case the 
probability that an arriving antineutrino mass eigenstate be detected as an electron 
antineutrino after crossing the Earth interior could be written as [73] : 

Pig ~ cos2^i3(l-p^), 

p2e - cos^^iapgj, (3.29) 

where is the two-flavor oscillation probability after crossing the Earth, can 
be computed under simple assumptions about the Earth interior. For example, if 
neutrinos travel only trough the Earth mantle and we assume a constant density, 
the oscillation probability is given by |75j : 

= sin' + sin' 29^, sin{2e^ - 29,,) sin' ("^^^^) , (3.30) 
where the in-matter mixing angle 6„i satisfies 13.131 

Figure 13.41 depicts the oscillation probability for antineutrinos as a function of 
the neutrino energy assuming different distances d traveled inside the Earth mantle. 
We have assumed standard values for ^ 4.5 g/cm^ and Ye = 0.5, and the 

mixing parameter values of sect. 13.1.21 

A more complex problem is to describe neutrino oscillations when neutrinos tra- 
verse the mantle and the Earth core. In a dedicated study Akhmedov jH5 obtained 
analytical relations assuming two layers of constant density. In that case the oscil- 
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lation probability p© for a neutrino which travel distances cIm + dc + dM could be 
parameterized in the following way [75] : 

= sin^ Ou + 0\ cos 2^12 + O1O3 sin 2^12 , (3.31) 
where the "oscillating" coefficients 0\ and O3 are the components of the vector 

= 2sin0jvf^^+ sin^ctc (3.32) 

having defined Y and as 



Y = COS0A/ COS0C — (Uf ■ ^c) sin0M sin^c* , (3.33) 
r ^ (sin2^^,0, -cos2^^), (3.34) 
_ sin 2^12 Am^2'^ 



sin 2^™ AE 



(3.35) 



The subindices M and C for and t in ()3.32|) and ()3.33|) refer to the mantle and 
core values of 6m and d. 

The distances traveled by neutrinos inside the mantle and core are given in terms 
of the Earth radius i?© = 6371 km, the core radius tq = 3486 km and the zenith 
angle z: 



M = R [ —COS z — J r'^/R'^ — sin z ] (3.36) 



dc = 2R^rl/Rl-sm^ z (3.37) 
In figure 13.51 we show the probability p© for different energies and zenith angles. 



3.2 The neutrino signal formalism 

The theoretical relation between neutrino emission in a supernova and the detected 
signal involves three elements: i) the neutrino emission, ii) neutrino oscillations and 
iii) the detection process. 

With the conversion probabilities we can proceed to determine the way neutrino 
oscillations mix the fluxes and spectra of the different neutrino flavors. 
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At the source, the total flux of neutrinos a, Sa{E, t) can be written in terms of the 
luminosity La(t), the neutrino average energy Ea{t) and the spectral distribution 
Fl^{E-t) as 

S^{E,t) = ^FT\E-t). (3.38) 

At the detector these fluxes are mixed to construct the observed antineutrino flux 
according to the conversion probabilities ()3.28p : 



L-'Sl''\E,t)=P-,S^{E,t) + {P-, + P^)S^, (3.39) 

where L is the supernova distance and we have used the fact that for all purposes 
the and fluxes are identical Sx = S-p ^ S- (see sect. 12. 

Using the unitarity property of the conversion probabilities J2a = 1 we can 
express the detected flux ()3.39|) in terms of just one probability, namely the electron 
antineutrino survival probability p = Pyj^^ : 



L^S-f *(^,t) = pS-e + (1 -p) S^. (3.40) 

From the detected flux of electron antineutrinos we can straightforwardly compute 
the rate of events at a given neutrino detector: 



dEdt 



Nt J dE' *(£;', t) a{E') e{E') n{E, E') , (3.41) 



where (j{E) is the detection cross-section, Nt is the number of target particles in 
the fiducial volume, e{E) is the detection efficiency and TZ{E, E') is the energy 
resolution function that accounts for the uncertainties in the measurement of the 
neutrino energies. 



3.2.1 Neutrino Detection 



The statistical procedure that will be introduced in the next chapter is based on 
two assumptions about the neutrino signal: it must have a large number of events 
(several thousands) and both the energy and time measurements must be available. 
Both conditions are fulfilled for electron antineutrino events in water Cerenkov de- 
tectors and in scintillator detectors. This is the main reason why we concentrate 
our attention on z/g. 
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Electron antineutrinos are detected in water and scintillator detectors through the 
inverse /5 process: 

V,+p^n + e+ , (Ercact ^ l.SMeV). (3.42) 

An accurate parameterization of the cross-section for this process has been recently 
provided by Strumia and Vissanij 



O^uj'^eP e n) _ rp 77'-0.07056+0.02018 In -0.001953 In^ 



10-43cm2 



pE ^rO.07056+0.02018 In iiii -0.001953 In" ii'i; (3 43) 



where E^. = Ep — Qnp with Qnp = rrin — mp ^ 1.293 MeV and all the energies are 
in MeV. This parameterization is accurate on a wide range of neutrino energies. In 
figure IT^ we compare this cross-section with the approximation ap{ueP e^n) ~ 
9.52 X 10~^^ cm^ PeEe and the corresponding time-integrated spectrum. We see that 
using an accurate cross-section the neutrino flux is reduced as much as 20%. 

In heavy water (D2O), z/g can be detected through the charged current reaction off 
deuterium: 



Ue + d^n + n + e+ (E.eact = 4:.03MeV) (3.44) 

SNO is the only operational heavy water detector and therefore, with fiducial vol- 
ume much smaller than Super Kamiokande, the number of neutrino events detected 
through this process will not be enough large to apply efficiently our method. 

Other reactions (as described in details in ref. |HS|), including Ve absorption in 
oxygen (water detectors) and carbon (scintillator detectors), will produce typical 
number of events below 10% of the statistics from 77g j3 process. Additionally, some 
of these reactions have large energy thresholds or cannot provide energy and/or time 
informations, a central condition for the application of any method. 

In our simulations we have assumed 100% detection efficiency above the energy 
threshold for all the detectors. This idealized situation will not change much the 
signal characteristics. 

The uncertainty in the energy measurement is an important aspect in our simulations 
and their analysis. We have approximated the resolution function TZ{E, E') with a 
Gaussian distribution with mean E equal to the measured energy, and as standard 
deviation AE equal to the energy resolution, that can be parameterized as|66j: 

A.E I "e E 
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Figure 3.1: Contours of equal neutrino jumping probability across a resonance region 
with different mixing parameters {/\m?,6). Dotted lines correspond to a 50 Mev neu- 
trino while dashed lines are for 5 MeV. Contours at the left range from Pj = 0.02 to 
Pj = 0.1 and at right from Pj = 0.9 to Pj = 0.98. It is possible to recognize three 
regions where jumping probabilities have different values: one were adiabatic transitions 
(Pj ~ OJ take place (upper-right) , a second one for which there are non-adiabatic transi- 
tions (Pj ~ Ij (lower-left) and an intermediate one where Pj has intermediate values for 
almost all neutrino energies. Jump probabilities are computed assuming the simple density 
profile p{r)Ye = 2 x 10^ g/cm^(r/10^ cm)^ |^ 
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Figure 3.3: Schematic illustration of how transition probabilities are computed using 
jumping probabilities and following the evolution of neutrino state along the resonance 
regions. The example illustrated is for Vf. transitions and normal hierarchy. 
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Figure 3.4: Oscillation probability in the Earth mantle for different distances dM 
traveled. We used Amj^ = 8 x 10"^, 9^2 = 30°, = 4.5 g/crn^ and = 0.5. 
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Figure 3.6: Upper panel: comparison between the approximate (dashed line) and the 
accurate cross-section \3.4'A )- Lower panel: time integrated spectrum obtained when the 
approximate (dashed) and the accurate cross-section are used. 
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Detector 


Eth {as, ^e) 


Fiducial 


jydct 






mass 






(MeV) 


kton (D 


= 10 kpc) 



Cerenkov 


SK[0] 
(H2O) 
SNO[0] 


5 
4 


(0.47,0) 
(0.35,0) 


32 


5,900 


- 9,990 




H2O 






1.4 


260 


- 440 




D2O 






1.0 


80 - 


■ 160 


Scintillator 


KamLAND [0] 
(N12+PC+PP0) 


2.6 


(0,0.075) 


1.0 


240 


- 400 


Cerenkov 


HK[0] 
(H2O) 


5 


(0.5,0) 


540 


100,000 


- 170,000 




UNO [0] 


5 


(0.5,0) 


650 


120,000 


- 203,000 




(H2O) 












Scintillator 


LENA [0] 
(PXE) 


2.6 


(0.1,0) 


30 


7,500 - 


- 12,600 



Table 3.2: The relevant Vg detection parameters for present and proposed large volume 
detectors. In the last column we give a range for the total number of expected Vf, events 
from a Galactic SN at 10 kpc. The larger and smaller numbers correspond respectively to 
a supernova model with large spectral differences between 17 ^ and (Supernova model I 
in chapter \^ and of a model where smaller spectral differences are obtained (Supernova 
model II). As regards to neutrino oscillations, we have assumed NH and sin^ ^12 = 0.26 
which give a mixing I^g : U^, w 4 : 3. Only charged current reactions, that provide good 
energy and time informations, have been considered. 

The value of the adimensional coefficients qe and Be and other characteristics of 
several operational detectors as well as of a few planned detectors are presented in 
table EH 

3.3 Characteristics expected for the signal 

Starting from the detected neutrino rate ()3.50|) we can compute the properties of 
the detected signal. The most important informations are: i) the number of events, 
ii) the time-integrated spectrum and iii) the energy-integrated time-profile. 
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Normally the analysis performed on a neutrino signal focus on one of these prop- 
erties. In this section we will evaluate these quantities under different conditions, 
in the attempt of understanding the interplay between the different effects for the 
determination of the signal characteristics. 

3.3.1 Number of events 

The number of events expected in a given energy and time interval is computed 
integrating the total rate in the desired region: 



The counting of neutrino events provides very important global information on the 
signal. The total number of events, for example, is directly related to the total 
energy released in the supernova. In table 13.21 we present the total number of 
events for a future Galactic Supernova in present and future detectors. 

In our analysis the number of neutrino events in certain energy and time intervals 
are clue properties of a signal. For example, we are particularly interested in low 
energy events because they are very sensitive to a neutrino mass, and how many 
neutrinos of this type could be present in a signal becomes a very important piece 
of information. 

In figure 13.71 we present the distribution of neutrinos inside several interesting E-t 
regions. To compute these numbers we used two very different supernova emission 
models and two extreme cases for the neutrino oscillation patterns. 

We can see that almost in all cases the majority of events will have intermediate 
energies and will arrive at late times (t > 1). Only few percent of the neutrinos will 
have energies below 10 MeV, this means that in a signal with 10^ events, hundreds 
of neutrinos of low energy will be available, enough to provide a good sensitivity 
to the neutrino mass. This suggests that it may be appropriate for the purposes of 
mass related analysis to have energy thresholds below the 10 MeV level. Just few 
neutrinos with energies smaller than 5 MeV are observed in almost all cases implying 
that an improvement in the energy threshold below this level will not change too 
much the potential to neutrino mass measurements. 



CHAPTER 3. SIMULATION OF A SUPERNOVA NEUTRINO SIGNAL 




(3.46) 



Time-integrated spectrum 



46 



100 



50 



10 



2.5 



0.1 1 

t(sec) 



30-5% 
(12-4%) 



10 



. 6-3% . 


16-7% 


■.■: . :48-8;3%: V^ , r 


(6-5%) 


(18-13%) 


.' ■ (62-7-6'%)) ' 






0.4-2% 






(1-2.5%) 






0.01-0.05% 






(0.02-0.07%) 

1 1 



o 

o 

D 

1/3 



Figure 3.7: Distribution of events in a supernova signal. Fractions are computed using 
\3.4()^ - The ranges refer to results obtained with two different supernova emission mod- 
els. The first figure correspond to Supernova model I (see chapter \^ with large spectral 
differences; the second one to Supernova model II with small spectral differences. Two 
different extreme mixing schemes were considered: a moderate mixing obtained when NH 
is assumed (numbers in parenthesis) , and a complete spectral swap that can be produced 
with IH and sw? > 10^^. 



3.3.2 Time-integrated spectrum 

One of the most interesting properties of a supernova signal is tlie distribution of 
neutrino energies. The energy spectrum of neutrinos could bring us clues about the 
neutrino emission process and will allow us to study the effects that oscillations in 
the supernova mantle and Earth interior will have on the signal. 

The time-integrated energy spectrum computed from the total detected rate (j3.41|) 
is given by: 
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Substituting the detected flux 5^°* <^^^ in <fn/dtdE we can rewrite fE in tlie 
following way 



where F-^[E) and Fx{E) are the independent time-integrated energy fluxes of Ve 
and Vr^ respectively. In the last equation we can see that neutrino mixing have a 
"modulation" effect on the neutrino spectrum distortion. 

In flgure 13.81 we depict for two different emission models, and assuming differ- 
ent neutrino oscillations schemes. There are two interesting facts in the spectral 
distortion observed in the flgures. The flrst one is that in the absence of Earth 
oscillations the distorting effects introduced in the case of normal hierarchy can be 
mimicked with an inverted hierarchy, if the mixing parameter sin^ ^13 has a value 
in the intermediate region of flg. 13.11 The other observation is that Earth matter 
effects will produce only a mild distortion. 

3.3.3 Energy-integrated time profile 

The energy-integrated time proflle is computed from the detected rate (j3.41|) as: 



ft provides global informations on how the neutrino flux changes in time, regardless 
of the neutrino energies. 

Figure 13.91 depicts typical time-proflles obtained for different oscillation schemes. 

We notice that besides a global modulation effect, the oscillatory distortion of the 
flux that was clearly observed in the spectrum can not be discerned in the time 
proflle. In particular. Earth matter effects become almost undetectable. This is 
easy to understand since being indeed small. Earth effects are averaged out when 
the integration in the energy is performed. 

3.4 Generation of full statistics signals 

Many of the techniques that have been considered to analyze a signal from a future 
Galactic supernova use one or more of the properties described in the previous 
section. To obtain the number of events or the integrated signal proflles it is not 



pF^{E) + {l-p)F^{E) 
F^{E)-il-pm{E)-F^{E)], 



(3.48) 




(3.49) 
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Figure 3.8: Time integrated spectrum for two different supernova model and oscillations 
schemes. Upper panel: supernova model I: (1) NH, (2) IH, sin^ ^13 = 10~^, (3) IH, 
sin^ ^13 = 10~^, Earth matter effects, zenith angle z = 100 (path just into the mantle) (4) 
same as (3) hut with z = 180° (mantle+core+mantle). Lower panel: supernova model II: 
(1) NH, (2) IH, sin^ 6*13 = 10~^, Earth matter effects, z = 100 (path just into the mantle) 
(3) same as (2) hut with z = 180° (mantle+core+mantle) 
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necessary to generate a complete signal, and then testing these techniques is rather 
straightforward. On the other hand, if we want to perform a statistical analysis of 
the whole statistics in the signal, as is done in only a few of the proposed techniques, 
it becomes necessary to generate a complete realization of the signal in all its details. 

There are two different ways to generate a detailed synthetic supernova neutrino 
signal. The first one is an intuitive procedure where the various physical processes 
are simulated. Neutrinos of different flavors are generated using the emitted fluxes 
Sa{E,t). Then the "seed" sample is filtered first according to the conversion prob- 
abilities, where neutrinos change flavor depending on their energy, and after with 
the cross section and efficiency of detection. The final signal is obtained joining all 
the resulting events. 

Despite the apparent naturalness, this method has several defects. Firstly, although 
sampling the emitted fluxes seems to be rather simple, because they are not affected 
by the distorting effects of oscillations, designing reliable filters using the conversion 
probabilities, particularly in presence of an energy dependence or a strong oscillatory 
behavior, is a non trivial matter. On the other hand since the seed sample to which 
the filters are applied has a finite number of neutrinos, statistical fluctuations could 
be amplified by filtering, producing artificial enhancements or depletions of given 
parts of the signal. 

The second method consists in sampling directly the detected rate (Pn/dtdE ()3.41|) : 

^^^^M = Nr StiE, t) a{E) e{E). (3.50) 

The signal is constructed generating pairs of energies and arrival times {E, t) 
(with as given by eq. ()3.46|) ). If a non zero neutrino mass is assumed, the time 
coordinate of each event is shifted according to the time-of-flight delay ()1.4|) . To 
account for measurement uncertainties (which affect mainly the measured energies) 
the energy of each event is reshuffled according to a Gaussian distribution with mean 
equal to the generated value E, and dispersion equal to the resolution of the detector 

The synthetic signals used in the MC analysis performed to test our method were 
generated through this procedure. 

Although the detected rate 1)3.411) could be in principle a complex multivariate func- 
tion (see fig. I3.10p . this method is free of the statistical flaws that could arise when 
the flrst procedure is used. 

In figure ITTUl we have summarized all the characteristics of a synthetic signal. We 
have plotted the surface defined by the total rate in the most general case when 
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Earth matter effects are present, in the E-t plane we show a scatter plot of the 
events generated with such rate and on the E-rate and t-rate planes we show the 
time-integrated spectrum and the energy-integrated time profile. 




Figure 3.10: Neutrino detected rate (surface), time-integrated energy spectrum (curve on 
E-rate plane), energy-integrated flux (plane t-rate) and scatter plot (dots on plane E-t) for 
a signal generated with Supernova model I assuming IH, siv? 6is = 10~^ and z = 180° 
(mantle- core-mantle) . The profiles have been normalized properly to fit into the plot area. 



Other details of supernova signal synthesis are presented in Appendix ^ where we 
describe the computer tool that we have developed to generate and analyze super- 
nova neutrino signals. 
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Mass limits with Supernova 
Neutrinos 



In this chapter we will describe a new statistical method to constraint neutrino 
masses using a high statistics supernova neutrino signal. The method was described 
by the first time in a intuitive form in ref. |83j . 

We have organized this chapter as follows. The first two sections are devoted to 
present the main idea underlying almost all methods to measure or constraint neu- 
trino mass using supernova neutrinos. Then we describe synthetically several of the 
recently proposed methods including those that were used to constraint the neu- 
trino mass using the SN1987A signal. An outline of the new method and the main 
assumptions behind it is presented in section 14.31 Section 14.41 presents a rigorous 
description of the method and several of its mathematical and numerical proper- 
ties. Details about the numerical computation of the statistics and other quantities 
related to the method are presented in sections 14.51 and 14.61 

4.1 Time-of-flight delay of supernova neutrinos 

It was realized long time ago that valuable informations on the neutrino masses could 
be provided by the detection of neutrinos from a Supernova explosion j2Zl EH 1231101 • 

The basic idea relies on the estimation or measurement of the time-of-flight (TOF) 
delay Attof that a relativistic neutrino of mass and energy Ei, traveling a distance 
L would suffer with respect to a massless particle: 
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L V 

E \ f im?\ 

f L \ /10MeV\ / \2 ^^^^ 
At„, - 5.1m.ec . (4.1) 

Figure 14.11 shows typical values of the TOF delay and its dependence on neutrino 
mass and energy. 

A natural question arises: what could be a "benchmark" in the signal that could 
play the role of a "massless particle" to estimate the neutrino delays? Alternative 
answers to this question are at the root of different approaches used to exploit this 
basic idea. We will summarize several of them in the next section. 

4.2 Neutrino mass and the TOF delay 

In order to measure or constrain a neutrino mass using the TOF delay is necessary to 
identify some kind of "timing" information in the signal. Many different ideas have 
been proposed since the seminal work by Zatsepin |2Zj. Some of them have been 
even already used to obtain information on neutrino mass from SN1987A signal O 

nsmniiiziiiH!. 

We could classify all the methods in three general classes: 

• Methods that use the change in the time spread of the signal due to delays 
between the arrival time of high and low energy neutrinos. 

• Methods that rely on the occurrence of specific timing events in coincidence 
with neutrino emission. 

• Methods that make use of a detailed statistical description of the signal. 
4.2.1 Time spread of the signal 

A supernova neutrino signal has a relatively short duration. Most of the neutrinos 
are emitted in just few seconds and although the detailed evolution of the neutrino 
flux depends on the supernova model, the total duration of the signal, which is 
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Figure 4.1: TOF delay as a function of the neutrino mass and the supernova distance 
(upper panel) for different neutrino energies. Lower panel depicts the continuous depen- 
dence on the energy of the TOF delay for different values of the neutrino mass. 
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determined by the time scale of neutrino diffusion in tlie supernova core, can be 
predicted from basic physical principles. 

Neutrinos of different energies emitted simultaneously from the supernova core will 
arrive at different times. This effect will increase the time spread of the signal in 
an amount directly related to the value of neutrino mass and supernova distance 
according to ()4.H) . Therefore, knowing the expected duration of the supernova 
neutrino burst and measuring the observed time spread it is possible to constraint 
or measure the neutrino mass. 

Using this kind of approach the total duration of the signal from SN1987A was used 
to estimate model independent limits on the electron neutrino mass in the range 
20 — 30 eV |13]. More stringent limits were obtained using specific assumptions 
about the signal time structure jlHl UHl EZl HH] • 

More recently a method to measure the time spread of neutrinos emitted during the 
shock breakout in a future Galactic supernova has been proposed [221 • The short 
duration emission of neutrinos in this phase, determined by the time scale of the 
shock wave emergence in the supernova core, could allow to set limits at the level 
of 3 eV without make any additional assumption about the time structure of the 
neutrino emission. 

4.2.2 Timing events 

Another way to estimate a neutrino mass with a supernova signal is based on "tim- 
ing" events that take place simultaneously with the supernova explosion, and used 
as benchmarks to measure the neutrino delays. 

Gravitational waves are expected to be emitted almost in coincidence with the neu- 
trino burst [Sni E]- If the supernova is close enough to allow a clear detection 
of gravitational waves, the peak in the emission of this radiation could be used to 
measure the delayed arrival of neutrinos in the early phases of the explosion, and 
constrain neutrino masses at the level of 1 eV j3T]. 

Another possibility is the early formation of a Black hole in the supernova core that 
will abruptly truncate the neutrino emission. With a non zero neutrino mass the 
cutoff in the observed flux will not be turn off sharply, and a measurement of the 
time spread of the neutrino signal could allow to put limits on neutrino mass at the 
level of 1.8 eV ^\^. 
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4.2.3 Detailed analysis of a signal 

A non-zero neutrino mass will produce a global distortion in a supernova signal that 
would be very difficult to be mimicked by astropliysical effects at the source. Using a 
detailed model of neutrino emission and performing a complete statistical analysis, 
values of the neutrino mass could be constrained by subtracting the effect of the 
masses and measuring the agreement between the emission model and the signal. 

A very accurate statistical analysis performed under these lines of reasoning was 
done recently by Loredo and Lamb jlHl- They used a Maximum Likelihood anal- 
ysis to fit several types of neutrino emission models with the SN1987A data, and 
Bayesian reasoning to account for prior information about those models and to per- 
form statistical comparisons between them. They obtained the limit m-p^ < 5.7eV 
which is more stringent than any other limit coming from the analysis of SN1987A 
signal, and close to laboratory upper bounds on the electron neutrino mass. 

The method proposed here belongs to this kind of approaches. In the following 
sections we present a detailed description of the method and its properties. 

In table 14.11 we summarize the neutrino mass limits from a supernova signal, for 
SN1987A and for a future Galactic supernova. 

4.3 Basic description of method 

The method proposed has three basic characteristics: 

1. It is based on a neutrino- by-neutrino kind of analysis and requires the deter- 
mination of the time and energy of each event in the signal^ 

2. The method uses the full statistics of the signal, i.e. every neutrino event in 
the signal is used to construct the statistical estimator. 

3. It can be applied independently of particular astrophysical assumptions about 
the characteristics of the neutrino emission (time evolution of the neutrino lu- 
minosity and spectral parameters) and does not rely on additional benchmarks 
events for timing the neutrinos TOF delays. 

The method relies on two basic assumptions: 

^j7e events on scintillator and Cerenkov detectors are particularly well suited for this kind of 
analysis. Additionally, since Fe will provide the largest number of events on present and future 
detectors, we will focus on these events in the description of our method and in the tests performed 
to measure its sensitivity. 
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Table 4.1: Limits on the neutrino mass from supernova neutrinos, SN1987A and from 
a future Galactic supernova neutrino signal. 

Limit Description References 

SN1987A 

"T'Fe ^ 20 — 30 eV Model independent mass limits using the total 0^ 

duration of the signal. 
"^Fe ^10 — 26 eV Model dependent mass limits with assumptions [13 IS1 1^ HHI 

on the time structure of signal. 
niy^ < 5.7 eV Complete statistical analysis using different [IH] 

neutrino emission models and Bayesian 

techniques. 

Potential 



'TT'Fe ^ 3 eV Model independent limit using time spread 

after shock breakout 
nij;^ < 1.0 eV Gravitational waves as benchmark to measure j51j 

neutrino delays in the early phase of the 

emission 

mj!^ < 1.8 eV Time spread of the residual neutrino emission [541 153j 
after the abrupt cutoff of the flux by the 
formation of a Black hole 



1. The first and most important one is that inside the collapsing core neutrinos 
are kept in thermal equilibrium, by means of continuous interactions with the 
surrounding medium, and therefore are emitted with a quasi-thermal spec- 
trum. As was explained in chapter |21 this is a solid prediction of almost all 
supernova models and simulations and has also been confirmed by the dura- 
tion of the SN1987A signal that constitutes an evidence for efficient neutrino 
trapping within the high density core. 

2. The second hypothesis is that the time scale for the variation of the charac- 
teristics of the neutrino spectrum is much larger than the time lags induced 
by a non- vanishing mass (say, much larger than 5 msec, see ()4.1|) ). Also this 
assumption is quite reasonable, since it is a robust prediction of all SN simu- 
lations 1^ inni IHH 0] that sizeable changes in the spectral parameters occur 
on a time scale much larger than 5 msec (see fig. 12. 8p . 

According to the first assumption, a high statistics neutrino signal can be considered 
as a 'self timing' quantity, since the high energy part of the signal, that suffers only 
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negligible delays, could determine with a good approximation the characteristics of 
the low energy tail, where the mass induced lags are much larger. The second as- 
sumption implies that the time evolution of the spectral parameters as inferred from 
the detected sample will reproduce with a good approximation the time evolution 
of the neutrino spectrum at the source. No additional timing events are needed, 
and each neutrino, according to its specific energy, provides a piece of information 
partly for fixing the correct timing and partly for measuring the time delays. 

Figure 14.21 illustrates schematically the basic strategy of the method. 
4.4 Formalism for the signal analysis 

In real time detectors, supernova electron antineutrinos are revealed through the 
positrons they produce via charged current interactions, that provides good energy 
informations as well. Each z/g event corresponds to a pair of energy and time mea- 
surements {Ei,ti) together with their associated errors. In order to extract the 
maximum of information from a high statistics SN neutrino signal, all the neutrino 
events have to be used in constructing a suitable statistical distribution. The Like- 
lihood Function is probably the best option to perform the analysis of the whole 
statistics of a supernova signal. 

To compute a LF we need to specify a probability for each event in the signal 
and therefore we need to construct or select a suitable analytical model to describe 
the neutrino emission from supernova. However, we are interested in evaluating 
how much information about the neutrino mass is possible to extract from a signal 
irrespective of our detailed knowledge on the neutrino emission process. Therefore 
we will need to marginalize the extra information about the astrophysical description 
of the signal from the effect of a mass. We have constructed our statistical method 
on the basis of Bayesian principles of inference. In Appendix^ we present a simple 
introduction to Bayesian reasoning. The definitions, results and terminology used in 
this section to describe our method is presented there, and therefore a first reading 
of that Appendix is recommended. 

The power of Bayesian statistics to perform an analysis of supernova neutrino data 
has been demonstrated with the analysis of the SN1987A signal by Loredo and Lamb 
in ref. jlHj. Our method is similar in several aspects to the Loredo's analysis but the 
fact that we will analyze signals with thousands of events will require a somewhat 
different mathematical and numerical approach. 
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Figure 4.2: Schematic explanation of the method strategy. Using the observed neutrino 
energy distribution at high energy the best-fit spectrum at a given time is found (upper 
panel). The neutrino mass can be measured by shifting the events with a time delay 
according to the test mass and event energy, until the corrected flux reproduces what is 
expected from the extrapolation to low energies of the fitted spectra at each time (lower 
panel). 
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4.4.1 The Likelihood Function 

We identify a signal as a set of pairs of energy and time measurements V = {Ei,ti]. 
The Likehhood Function (LF) £, associated to a particular parametric model of 
the emission ^A = M.{{6}) and a given neutrino mass m^, could be schematically 
written as: 

p{V\m\ M) = C{m\ {9}- V) = \[ f{E,, t,) , (4.2) 

i 

where the index i runs over the entire set of events in the signal. Here, f{E,t) 
represents a probability distribution function (pdf) used to evaluate the contribution 
to the likelihood of a single event. 

Using the emission rate expression in eq. (j3.38p and the total detected rate in eq. 
(|3.41|) . / can be written as: 

fit, E- m\ {6}) = 4>{t) X F{E; t) x a{E) , (4.3) 

where the normalization constant is = j dt J dE (pit) F{E;t) <y{E). In the above 
definition we identify three components: the neutrino time profile flux (p{t), the 
energy spectrum F{E; t) which in general evolves in time, and the detection cross- 
section (j{E) which is a well known function of the neutrino energy. 

4.4.2 Posterior probabilities 

In Bayesian inference the likelihood of a given model in the presence of some ex- 
perimental evidence must be always weighted with the prior probability that we 
could give to the model itself and to the measurement process (see Appendix 
Therefore, when using Bayesian principles the LF is just the first step to obtain the 
posterior probability for the theoretical model. 

Using Bayes theorem ()A.4jl the posterior probability of a signal model is given by: 



2 ,n^^^^ _ C{m\ {e}-V)p{m'')pm) 



p{m\{e]\V) _ ^^^^ 
p{rr?,{e}\V) = C{m\{e};V)p{m^). (4.4) 

In the last equation the evidence, p{T>), which does not depend on the parameters 
of the model {6}, has been absorbed in a normalization constant A/". Flat prior 
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probabilities for the model parameters p{{0}) = const, (as it is often the case) can 
also be absorbed in the normalization. 

The choice of some suitable form for the prior probability for is more subtle. We 
have used for simplicity a step function p{m^) = G(m^) = 1, (0) for > 0, (< 0) 
to exclude unphysical values of the neutrino mass. It must be noticed however that 
this is not the better way to describe our present knowledge about neutrino mass, 
since a prior probability flat on implies an uneven prior probability on m. We 
will discuss later in sect. 15.41 how the choice of other priors can affect the results of 
the analysis. 

The posterior probability in eq. ()4.4|) corresponds to the joint probability of the 
emission model characterized by the set of parameters {6} and a given neutrino mass. 
However, here we are just interested on what the signal could tell us about the mass 
and then we need the posterior probability just for this quantity. Bayesian statistics 
provides a natural way to marginalize the nuisance parameters (the emission model 
parameters {6}) by integrating them out 



Our final goal will be to compute p{ml\V) and to derive from it all the possible 
probabilistic informations contained in the signal. 

All methods to constrain the neutrino mass using a supernova signal are directly 
sensitive to m^. However it is common to search for a way to express the limits on 
ml as limits on m^. In our case that translation is natural since knowing the pdf 
for ml, to write down the pdf for m^ requires only a variable transformation. The 
posterior pdf for mi, reads: 



4.4.3 Best-fit value and neutrino mass limits 

Once the posterior pdf for ml and are computed it is possible to derive proba- 
bilistic statements about the neutrino mass. There are three pieces of information 
which can be obtained. 

Firstly the LF alone could provide the best-fit values for ml and for the parameters 




(4.5) 



p{m\V) = 2|m|p(m^|P) 



(4.6) 
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of the emission model (the ones that maximize the likelihood). This Maximum Like- 
lihood (ML) analysis gives us the most probable value m^^ and its error, compatible 
with the neutrino signal, but conditioned to the specific set of best-fit values {6'fit} 
of the model parameters. 

For practical reasons in the ML analysis the logarithm of the likelihood function, 
the log-likelihood function (log-LF) is used. Using eq. ()4.2j) the log-LF reads 



log/: = J2\ogf{Ei,U;{9}) 

i 

= J2 log X F{Ei; U) X a{Ei)] - log AT. (4.7) 

i 

It is the maximum of the posterior pdf p(m^|D) which actually provide us with the 
most probable neutrino mass, irrespective of the value of the model parameters and 
taking into account all the prior information. 

The most useful information about the mass is obtained computing the credible 
regions (CR) that are the range where the neutrino mass lies with a certain proba- 
bility. 

For our purposes we will use three CR that respectively give us the lower and upper 
limits in m^, and an upperbound on niy. 

Upper and lower limits for m^. The upper and lower limit on is obtained 
requiring that 



p{ml\D) dml = CR, [U 

-oo 



and 



p{ml\D) dml = CR. (4.9) 

w 

Figure 14.31 illustrate the definition of these limits. 

Upper bound for m. If the neutrino mass is too small to be discerned an upper- 
bound on m can be obtained from its posterior pdf (|4.6|) requiring 

/•rriup 

/ p{m^\D)dm^ = CL, (4.10) 
Jo 
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Figure 4.3: Schematic illustration of the determination of mass limits from the posterior 
probabilities for rn^ (upper panel) and (lower panel). 
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4.5 Construction of the Likelihood Function 

To compute the LP ()4.2p two basic elements are required. First, we need to deduce 
the time evolution of the spectral parameters from the sample in order to describe 
the neutrino energy spectrum at any time. Second, we have to find a suitable class 
of parametric analytical function that could fit reasonably well the detected flux. 

4.5.1 Neutrino spectrum 

According to the first assumption in sect. 14.31 the spectrum can be reasonably de- 
scribed by a quasi-thermal (analytical) distribution. If for example a distorted 
Fermi-Dirac function ~ [1 + exp{E/T — ri)]~^ is used, as was done in ref. jH31; F 
can be parameterized in terms of a time dependent effective temperature T and de- 
generacy parameter rj describing the spectral distortions. According to the second 
assumption, the time dependence of the relevant spectral parameters can be inferred 
directly from the data. 

Here, we model the observed neutrino spectrum F{E] t) by means of an a-distribution 

(EH), 

F(E,e(t),a(t)) = N{e,a) (E/e)%-("+^)^/% (4.11) 

where e(t) and a{t) describe the evolution in time of the spectrum. 

The choice of this distribution to fit the observed neutrino spectrum instead of the 
physically better motivated Fermi-Dirac, obeys similar reasons than those presented 
in sect. 12.2.31 Starting from a discrete sample of neutrinos, a Fermi-Dirac spectrum 
can be reconstructed only by carrying out numerical fits to the energy momenta until 
the correct values of T and rj are determined through a minimization procedure. 
In contrast, the a-distribution can be straightforwardly determined through the 
simple analytical relations (j2.8|) connecting the spectral parameters and the two 
first momenta of the energy distribution. 

Since the effect of the detection cross-section (see sec. 13.2.11) modifies the measured 
energy distribution, we need to find a way to reconstruct the energy momenta of 
the emitted distribution from the set of observed energies. 

If we call Fe{E) the energy spectrum at the source (where for simplicity we have 
dropped the temporal dependence) and Fii{E) the (unknown) observed neutrino 
spectrum, Fii{E) = Fe{E) x (t{E)^ the n*/i-momentum of is given by: 
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j^dEF^{E) 
J^-dE E-F,{E)/a{E) 



j^dEF,{E)/a{E) ■ 
The last relation allows us to define the new cross-section weighted momenta: 

This weighted momenta provide us an estimate of the momenta of the emitted energy 
spectrum and allows us to find the spectral parameters around a given time. So, in 
order to obtain the spectral evolution, we slice the sample in many time windows 
and find the values of a and e inside them. 

A couple of words about the procedure of signal slicing procedure are necessary to 
explain in more details how the functions giving the spectral parameters evolution 
are constructed. 

Any shcing procedure faces the effects of statistical fluctuations. In our case these 
fluctuations come from the limited statistics on each window and from the fact that 
neutrinos from different times and with different spectral temperatures at the source 
could be mixed up inside a single window. 

To reduce statistical fluctuations in the spectral parameters estimation the construc- 
tion of the time windows and the processing of the resulting information has been 
performed in the following way: 

1. The width of each window is chosen large enough to contain a sufficient number 
of neutrinos (a few hundreds). This reduces the statistical uncertainties coming 
from finite statistics effects. 

2. Time windows are overlapped: the central value of each new window is de- 
termined as tn+i = tn + St, with St <S Attof- Therefore many neutrinos of 
neighbor windows will be contained on each window and averaging the result- 
ing fluctuations will produce smooth estimates for the spectral parameters (see 



Once the spectral parameters in each window have been estimated their val- 
ues are smoothed to wash-out the statistical fluctuations and the effect of 
neighbors. 
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4. The time variation of a and e is finally fitted with a non-linear least-square 
analysis. 

A schematic illustration of the above procedure is depicted in fig. 14.41 

In figure l^31 we compare the estimation of the spectral parameters with the described 
procedure, for a set of synthetic supernova signals and the real values of those 
parameters used to generate the samples. 

4.5.2 Detection cross-section 

Also for the analysis we have used the accurate parameterization presented in 
sec. 13.2.11 This is the same cross-section used to generate the signals. 

4.5.3 Flux model 

Coming back to the problem of constructing the Likelihood function, and in partic- 
ular of choosing a specific time profile for the neutrino flux (namely the model M) 
we have proceeded according to the following requirements: 

1. the analytical flux function must go to zero at the origin and at infinity 

2. it must contain at least two time scales: the initial, fast rising phase of shock- 
wave breakout, and the later Kelvin-Helmholtz cooling phase. Another pa- 
rameter could be required to describe the accretion phase, i.e. the transition 
point between these extreme phases. 

3. it must contain the minimum possible number of free parameters to avoid 
degenerate directions in parameter space and to speed up the Maximum Like- 
lihood analysis. Still, it must be sufficiently 'adaptive' to fit in a satisfactory 
way the numerical flux profiles resulting from different SN simulations, as 
well as flavor mixed profiles as would result from neutrino oscillations (see 
section 13. 

In order to test how much our results on the neutrino mass will depend on the 
specific flux profile, we have performed several tests using two different flux models 
constructed following the above criteria (results are described in sect. I5.2.2j) . 

In the following paragraphs we will describe some of the main properties of these 
two general fluxes. 
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Figure 4.4: Schematic illustration of the spectral parameter estimation using overlapped 
time windows. The signal (upper panel) is sliced in overlapped time windows with enough 
statistics (different points belong to three selected time windows). The emitted spectral 
mean energy e is estimated using the cross-section weighted momenta of the events con- 
tained on each window. The result is the polygonal line in the lower panel. The final values 
used as parameter e in the analysis are the smoothed values of the dotted line. 
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Figure 4.5: Results of the spectral parameters estimation for 40 synthetic signals. Upper 
panel (lower panel) : spectral mean energy e (spectral pinching a), estimated from the sam- 
ple with the procedure described in the text (oscillatory curves). Continuous line represent 
the expected value after the mixing ofVf, andv^ fluxes due to oscillations in the SN mantle. 
The acceptable level of agreement confirms the validity of the procedure. 
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Flux model I: exponential rising and power-law decay 

The following model for the flux, in spite of being very simple, has all the required 
behaviors, and moreover it showed a remarkable level of smoothness and stability 
with respect to numerical "extremization" and multi-parameter integrations: 

This model has five free parameters that on the l.h.s of ()4.13|) have been collectively 
denoted with {6}: a time scale ta for the initial exponentially fast rising phase, a 
second one tc for the power law cooling phase, two exponents Ua and Uc that control 
the detailed rates for these two phases, and one additional exponent Up that mainly 
determines the width of the "plateau" between the two phases (see figure . 

For the purposes of our numerical tests we have reduced the number of parameters 
to three by fixing the values of exponents Ua = 1 and Uc = S found with a fit of 
realistic time profiles. Although for each supernova model it is necessary to find the 
best values of these exponents, the election of these parameters will not affect too 
much the results of the analysis. 

Flux model II: Truncated accretion and power-law decay 

Here we describe a second analytical model for the evolution in time of the emitted 
neutrino flux that is based on the phenomenological analysis presented in sect. 12.1.41 

The analytical model in (|2.4j) does not fulfill all the requirements: for example, it 
does not include a description of the shock breakout phase and therefore, does not 
go to zero when t 0. To correct this we have included to ()2.4|) a multiplicative 
damping factor. 

Our phenomenologically motivated flux model reads: 



0(t; {9}) 



(t/U)''^ it- 

(tc/tr it- 



(1 + t/tfe)"" 

.0). 

> oo). 



C 



t/tc 



(4.14) 



This model has 11 parameters: four time scales, td for the initial fast rising, ta and 
tb for the flux plateau and the early flux decay and tc for the long term decay; four 
exponents, to which controls the rate of each component; and two amplitudes 
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Figure 4.6: Flux model I i4-l^ o.'i^d the effect of each parameter. Upper panel: the effects 
of parameters 1 (time scale tc) and 2 (rate of decay Hc); as can he seen the effect of both 
parameters are very different and little correlation between them is expected. Lower panel: 
the effects of parameters 3 (time scale ta) and 5 (rate of exponential rising Ua); also in 
this case little correlation is expected. 
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A and C which determine the relative contribution of the accretion and coohng 
terms to the total flux. 

As was done with the more simple flux model, we need to reduce the number of 
free parameters to simplify the extremization procedures. Studying the impact of 
each parameter (see flg. 14. 7p we have determined that enough flexibility could be 
still obtained free just three relevant parameters: the cooling exponent Uc, and the 
time scales of fast rising and early decay (truncated accretion) t^. The rest of 
parameters could be flxed at the beginning performing a simple flt of the signal time 
profiles with the flux model. 

In presence of neutrino oscillations the detected flux will be an admixture of and 
fluxes. Since the mechanism responsible for the emission of the two species of 
neutrinos are different, a simple description of the resulting flux is not possible. 

Given that in the Likelihood analysis we will set the origin of times in coincidence 
with the flrst neutrino detected, and this obviously cannot correspond to the origin 
of time of a flux function which must satisfy (j){0) = 0, an additional parameter 
6t is needed to allow the function to shift freely along the time axis according to 
cj){t) ^ cj){t + dt). 

4.5.4 LF regularization 

Until now we have not discussed the role of the mass in the likelihood function. 

The neutrino mass determines the value of the time shift At = At(m^, E, L) applied 
to the event probability ()4.3p in the deflnition of the LF ()4.2|) . In that sense we could 
think about the neutrino mass as another parameter of the pdf /. 

However, in order to make more transparent the way the neutrino mass enters into 
the likelihood, we could treat the effect of this parameters in the following completely 
equivalent way. For each new value of the tested was we could think as if the sample 
were changed: all neutrinos are shifted to smaller times according to At{ml, E, L) 
and we deflne a new set of pairs of energies and times {Ei,tf^) (see flg. 14. 8p . The 
LF of the emission model for this particular mass will be, 

C = l[f{E,,tt;{9}). (4.15) 

According this procedure, it is natural that some neutrinos can end up with negative 
values tf^ < 0. This effect would induce a compensation through a time shift of the 
flux model. The procedure is illustrated schematically in flgure 14.81 
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Figure 4.7: Behavior of the flux model II in eq. i4.15{ ) under changes in its parameters. 
We start with the reference set of parameter values: (1) A = 0.2, (2) C = 1.5, (3) ta = 0.2, 
(4)tb = 300, (5) tc = 0.1, (6) td = 0.05, (7) Ua = 0.3, (8) = 0.4, (9) = 2.0 and (10) 
Ud = 1.5. In each panel a pair of parameters are studied. As can he seen in panels \J7/jf i,b,c 
the parameters A, C , ta, Ua and tc have similar overall effect on the flux shape. To avoid 
correlations we have left free just one of them (ta). On the other hand Uc, td and Ud 
determine different behaviors and therefore are all left free. 
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Figure 4.8: Schematic illustration of the two-step process used to evaluate the likeli- 
hood function of a given model. Given a test neutrino mass a new sample (asterisks) is 
constructed shifting the times of all the events. The LF is computed by calculating the 
probability f of the events in the shifted sample, by applying a global shift to the flux model 
and spectral parameters (curves). 
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Because of the finite energy resolution of any detector, the measured neutrino en- 
ergies that are used to evaluate the time shifts for neutrinos in the sample do not 
correspond to the true energies that determine the real neutrino delays. Therefore, 
even when the correct value of the test mass is used, the time-shifted neutrino sample 
will not correspond exactly to the sample originally emitted. Although completely 
natural (as well as unavoidable), this behavior can produce a dangerous situation. 
When the energy measurement yields a value smaller than the true energy, a neu- 
trino arrival time can be shifted to negatives values where the flux function, even 
with a non-zero time shift, vanishes, implying that the log-Likelihood diverges. This 
would imply rejecting the particular neutrino mass value for which the divergence 
is produced, regardless of the fact that it could actually be close to the true value. 
Figure illustrates this point. 

To correct this problem we have adopted the following procedure. The contribution 
/ to the Likelihood ()4.15|) of a neutrino event with measured energy Ei ± AEi 
for which, after subtracting the delay 5ti = mlL/2Ef , we obtain a negative value 
t^'^ < (or a value close to the origin of the flux function ti ~ —6t) is computed by 
convolving it with a Gaussian Q(t; tf", cr*) centered in tf^ and with standard deviation 
a, = 2 5UAEi/Ei : 



Clearly this regularization of the divergent contributions to the Likelihood is phys- 
ically motivated by the fact that the origin of the problem is the uncertainty in the 
energy measurements, that translates into an uncertainty in the precise location in 
time of the neutrino events after the energy- dependent shifts are applied. 

The regularization procedure is one of the novel features of the method proposed. 
Although it is not free of numerical problems (e.g. the likelihood can have dis- 
continuities that must be corrected as illustrated in fig. I4.1U|) . the most interesting 
property of this procedure is that of allowing the inclusion in the likelihood of infor- 
mations about the resolution of the detector: the same event in different neutrino 
detectors will have a different contribution to the LF, when shifted to the forbidden 
region t^^ < —6t), because the width of the resolution function Q in fl4.16|) will be 
also different. Therefore detectors with a better energy resolution will be more sen- 
sitive to changes in the position of the events and therefore in the tested neutrino 
mass. 




(4.16) 
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Figure 4.9: Schematic illustration of the problem introduced because the detected neutrino 
energies are not the true energies, after shift according to a test mass. The asterisks 
represent the detected events. Triangles are at the original energy of the event. Bars 
indicate the size of the energy resolution of the detector (Super Kamiokande resolution 
is used). When shifted with a test mass the position of the detected events are different 
from the original positions. This is is a natural result of the energy dependence of this 
shifts. However, when the measured energy is smaller than the original one, as in the 
lowest energy event, even if a mass equal to the real neutrino mass is used, the resulting 
time shift could bring the event to a region where the log-LF diverges. 
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Figure 4.10: Time profile of f for a 10 MeV neutrino using three different schemes 
of regularization: no regularization at all (continuous line), regularization of events just 
into the forbidden region t^^ < —5t (triangles) and regularization of events inside and 
outside (but close) the forbidden region (dashed line). The last strategy intends to avoid 
discontinuous changes in the likelihood produced when a single neutrino event leaves the 
region where f is regularized. 
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4.6 Likelihood marginalization 

The LF alone cannot give us the information that we want about the neutrino mass 
without reference to the other parameters used to describe the signal. Indepen- 
dent probabilistic information on is obtained by calculating the posterior pdf 
p{ml\V) that results from marginalizing the LF with respect to the nuisance (flux) 
parameters. 

In practice, the marginalization procedure can be a very expensive numerical mul- 
tidimensional integration of a particularly "heavy" function. Each evaluation of 
the LF for a typical signal from a future Galactic supernova could require several 
thousands of evaluations of relatively complex functions (flux model, spectrum and 
cross-section). With all these conditions it is clear that the CPU time required to 
carry out all the necessary operations would be exceedingly large, specially consider- 
ing that to draw any statistical conclusion about the quality of our method, we need 
to analyze a large set of neutrino samples, corresponding to different SN models, 
SN-earth distances and also to different detectors. 

A way to avoid this problem, is to approximate the marginal posterior probability 
with the profile likelihood (PL) C{ml;D), that corresponds to the trajectory in 
parameter space along which for each given value of ml the Likelihood is maximized 
with respect to all the other parameters (see fig. I4.11|) . 

It can be shown that for a multivariate Gaussian the PL coincides with the marginal 
posterior p{ml\D) (for details see appendix , and therefore our results will be re- 
liable to the extent the Likelihood approximates well enough a normal distribution. 
In fig. 14.121 we compare different contours in parameter space for log£(m^, {6}; D) 
with those of a corresponding multivariate normal distribution, with the same mean 
and covariances than the likelihood. We see that within the region where the con- 
tributions to the integrations are large, the behavior of the Likelihood is indeed 
approximately Gaussian. 
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Figure 4.12: Contours of the log-LF compared with the contours of a Gaussian dis- 
tribution of the same mean and covariance in four different two parameters spaces: 
versus \J.l^ : the time shift 6t of the flux function; the time scale of the fast initial 
rase ta (n^ = 1 ); \4-l^ : the time scale of the cooling phase t^ (n^ = 8); \4-l^ : the ratio 
^c/np (see eg. \4-l^ - Por the construction of this contours we have used a particular 
parameterization of the flux-model \4-l'J\^ , the same used in ref. \8'J^ . in order to minimize 
correlations between parameters. The contours are plotted at 0.5, 1.0, 1.5 and 2.0 a. 
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We present in this chapter the results of a series of tests performed to study the 
behavior of the method and to estimate its sensitivity to a neutrino mass. 

The tests were performed with a wide range of emission (SN model), propagation 
(neutrino oscillations) and detection conditions, as described in sect. 15.11 The gen- 
eral results of the ML analysis are summarized in sect. 15.21 To quantify the sen- 
sitivity of the method we have measured two statistical properties for each set of 
analyses. The results of these measurements are presented and analyzed in sect. 15. 3| 
Finally, we discuss in sect. 15.51 the implications of the results and its relevance in 
the context of current and future neutrino mass searches. A summary and the 
conclusions of this work are presented in this section. 

Part of the material presented in this chapter were published in a similar form in 
the paper [Sn] 

5.1 Description of tests 

In the absence of a real high statistics supernova signal, we can only test our method 
applying it to realistic synthetic signals generated with the methods described in 
chapter |21 However of course we cannot derive any general and reliable conclusion 
from the analysis of a single signal and just one particular emission model and 
supernova set of properties. 

In order to have an overall understanding of what to expect if our method could 
be applied to a real signal, we have carried out a Monte Carlo (MC) analysis using 
many synthetic signals, that were generated under a wide range of conditions for the 
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emission (SN fluxes and spectra, SN distance), propagation (oscillations and TOF 
delays) and neutrino detection (number of events, energy resolution and energy 
thresholds) . 

We divide our MC analysis in a set of different tests. For each test we applied the 
method to a sample of ~ 40 synthetic signals each one generated with the same set 
of input conditions (we can imagine the samples conforming one test, as a set of 
many realizations of the same supernova, detected with the same detector). 

The input conditions considered to construct the whole set of tests were: 

• Neutrino emission. Two different supernova models were used: 

the first model [supernova model T) corresponds to the results of a simulation 
of the core collapse of a 20 Mq star which was performed with the code devel- 
oped by the Lawrence Livermore Laboratory group jHEl- These results were 
used to study supernova r-processes [Hni- Figure I^TTk depicts the neutrino lu- 
minosity and mean energy evolution obtained in this simulation. This model is 
characterized by large spectral differences between electron and non-electron 
neutrino species, which, as explained in sect. 12.31 is a result of the approximate 
description of the transport in the SN core. 

Neutrino spectra and their evolution were not published in the original paper 
by Woosley et al. jBU]. We have therefore used the spectral shapes taken 
from the dedicated study of Janka and Hillebrandt jHZ| and we have assumed 
for simplicity that the a-parameter, which quantify the spectral distortions, 
remains constant during the emission. We have used the standard values 
a-p^ = 3.5 and a^^ = 3.0 for this quantity. 

The second model {supernova model IT) correspond to the recent state-of- 
the-art hydrodynamical simulation of the core collapse of 15 Mq progenitor 
star 121 EI] carried out by means of the Garching group code [HHlEniinni- This 
simulation included a more complete treatment of neutrino opacities PTIIH^ IH]. 
that resulted in smaller spectral differences between 17^, Ve and (see fig. 15. lb ). 

The Garching group simulations were stopped after 750 msec, and the results 
were not completely reliable already after the firsts 300 msec jHl]- Since our 
study assumes a signal time duration of 20 sec, after which the neutrino flux is 
assumed to become undetectable, we have done a conservative extrapolation 
of the signal to later times. For the luminosities we have assumed a power 
law decay in agreement with general results of SN simulations |56| IHH IHOl 12] 
while for the mean energies we have assumed a mild decrease and a constant 
pinching after 750 msec. 
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(a) Supernova model I ( adapted from (b) Supernova model II ( adapted from 
ref. JEBl) ref. JES) 

Figure 5.1: Luminosity, mean energy and a-parameter for andVx in the supernova 
models used in this work to generate the synthetic signals. 
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We have also varied SN-earth distance into the representative values 5, 10, 
and 15 Kpc. 

• Neutrino oscillations. To perform the MC analyses we used an oscillation 
scheme with a normal mass hierarchy and the LMA mixing parameters set (see 
sect. I3.1.2|) . With this oscillation scheme the observed flux of Ve correspond 
to an admixture of ~74% and ~26% V^- 

The case of a complete swap ^ V^. that can result from an inverted mass 
hierarchy and |t/e3p ~ 10~^, can be considered equivalent to the non-oscillation 
case (studied in ref. jH21) with a larger average energy for 77g. Indeed we can 
consider the results of the analysis performed with supernova model I as the 
harder spectrum version of that corresponding to supernova model II. 

We do not include Earth matter effects, since they will depend on the specific 
position in the sky of the SN relative to the Earth, on the specific location of 
each detector and on the time of the day. However, given that even with a 
dedicated analysis it appears quite challenging to identify clearly these effects 
[7T] . we believe that this neglect is of no practical importance. 

• Detection. As explained in the previous chapter the application of our 
method requires the knowledge of both energy and time of each neutrino event. 
On the other hand the spectral estimation procedure described in sect. 14.51 
requires a large statistics of neutrinos. Electron antineutrino detection will 
provide by far the largest number of events from a future Galactic supernova, 
and therefore our analysis will be restricted to this kind of neutrinos and to 
detection processes capable to measure the energies and arrival times with 
good precision. 

Today, the Super-Kamiokande (SK) detector, having the largest fiducial vol- 
ume among all operative detectors, suited for the kind of analysis we propose. 
Operative scintillator detectors are characterized by a lower threshold and/or 
better energy resolution, which in principle represents an advantage in the 
analysis. However, their significatively smaller fiducial volume imply the pos- 
sibihty of detecting only a few hundred of events (see tab. 13. 2^ . Such a low 
statistics is not well suited for the application of our method. However, since 
KamLAND (KL) is located at the same site as SK (and equally affected by 
neutrino oscillations in the Earth) KL events can be the high statistics of SK, 
possibly yielding some sensitivity improvement. We have performed this kind 
of joint analysis to study the impact of the addition of KL events to the SK 
signal. 



CHAPTER 5. RESULTS AND DISCUSSION 



5.2. RESULTS OF THE MAXIMUM LIKELIHOOD ANALYSIS 84 



For the future large volume detectors we have chosen two of the most interest- 
ing proposals: Hyper-Kamiokande (a megaton water Cherenkov detector) 
and the multi-kiloton scintillator detector LENA (Low Energy Neutrino As- 
trophysics) jns]- 

The properties of all the detectors used in the MC analysis are summarized in 
tab. 13.21 in sect. 13.31 

The following table defines the set of labels that will be used hereafter to identify 
the different tests of the method and their particular input conditions: 



Label 


Supernova 


Detector(s) 


Distance 


Flux 




model 




(kpc) 


model 


TSNI-1 


model I 


SK 


10 


model I 


TSNI-2 








model II 


TSNI-3 






5 


model I 


TSNI-4 






15 




TSNI-5 




SK + KL 


10 




TSNI-6 




LENA 






TSNI-7 




HK 






TSNI-8 




HK 




model II 


TSNII-1 


model II 


SK 


10 


model II 


TSNII-2 








model I 


TSNII-3 






5 


model II 


TSNII-4 






15 




TSNII-5 




SK -1- KL 


10 




TSNII-6 




LENA 






TSNII-7 




HK 







In the last column of the table we have indicated what flux model (model I ()4.13|) . 
model II (j4.15|) was used to compute the signal pdf (j4.3|) . 

For the signals generated with supernova model II we have mainly used the more 
flexible flux model II. The analysis with such a flux in the case of supernova model 
I does not show any relevant difference with respect to the simple flux model 1. 

5.2 Results of the maximum likelihood analysis 

To illustrate some of the most interesting properties of the method and before pro- 
ceeding with the estimation of its sensitivity, we performed a set of MC analyses 
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using an hypothetical MC neutrino mass = 1 eV. 

Several interesting facts about the behavior of the method could be observed when 
the output of the ML analysis, namely the best-fit values of mass and flux parame- 
ters, are studied. In the following two sections we describe these results. 

5.2.1 Best-fit values and limits on ml 

The best-fit values for ml and its limits are summarized graphically in the band-plots 
presented in figs. 15. 2115. 51 

All the signals on a band plot were generated under the same set of input conditions. 
The squares indicates the position of the best fit value and the "error bars" 
correspond to the lower and upper limits computed from the posterior pdf p{ml\V) 
using (jOl) and (Oil . 

Some of the analysis have exceedingly large error bars when compared with the rest 
of analysis in the same set. In most of the cases they correspond to badly behaved 
analysis, where the extremization procedures used to compute the profile likelihood 
found peculiar directions in parameters space corresponding to large correlations 
between the mass and the flux parameters. In a few cases such "failures" can be 
caused by just one or two neutrinos that randomly get abnormal values of their 
energy. Exclusion of these neutrinos can generally rescue a well behaved analysis. 

The first remarkable fact is that regardless of the wide range of input conditions (su- 
pernova emission model, particular detector, flux model used to analyze the samples) 
in almost all cases the average value of ml is close to the MC mass used to generate 
the signals. This is a clear indication that the signal can act as a self timing ob- 
servable and that no particular assumption on the signal time structure is needed to 
measure the value of the mass. This fact represents a confirmation of our original 
hypothesis. 

The only feature of the analysis sensible to particular input conditions is the disper- 
sion around the best-fit values of the mass. As expected, a lower statistics increase 
the dispersion and reduces the sensitivity to a neutrino mass. This can partially ex- 
plain the slight difference between the results for the two supernova models. Since 
supernova model 1 signals have a larger number of events than supernova model II, 
the dispersion is smaller for the first model. We can conclude that a harder neutrino 
spectrum would increase slightly the sensitivity to a a neutrino mass, by increasing 
the statistics of the signal. 

Results obtained with just SK events (TSNI-,II-1) and SK-I-KL events (TSNI-,II-5) 
do not exhibit any appreciable difference, meaning that the better energy resolu- 
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Figure 5.2: First set of band plots for tests performed with supernova model I and 
"^MC ~ 1-0 eV^. The best-fit values for rn^ obtained from the ML analysis are represented 
by the filled squares. The error bars extend from inf^^ to m'^p. 
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Figure 5.3: Second sets of band plots for tests performed with supernova model I. 
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Figure 5.4: First set of band plots for tests performed with supernova model II. 
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Figure 5.5: Second set of band plots for tests performed with supernova model II. 
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tion and lower threshold of KamLAND cannot compete with the SK much larger 
statistics. 

A qualitative overall analysis of these results shows that sensitivities at the level of 
1 eV level can be achieved with SK and future LENA detector, while the sensitivity 
will be improved by a factor of 2 with a future megaton water Cerenkov detector. 



5.2.2 Fit to the flux model 



The ML analysis of the signals described in the previous section provides in each 
case also a set of best-fit values of the flux model parameters. The fit to the flux 
model is actually of little interest for the purposes of our work, but a comparison 
between the results of such a fit with the input time profiles and with the histograms 
of the detected samples, reveals other interesting aspects of our method. 

It must be realized that a direct comparison between the flux model and the 
observed time profile cannot be performed directly. The energy-dependent 

detection cross-section depletes the flux at early times when the mean energy of the 
neutrinos is lower. The comparison must be done between the observed profile and 
the energy integrated time-profile given by 



N- 



dE F{E-t)a{E) 



(5.1) 



In figures l^^5 . 71 we depict comparisons between ff°^{t) and fl^^^^^it) for several of 
the tests. Additionally we have included there the expected time-profile computed 
from the signal rate (j3.4H) : 



The last comparison allows us to understand how well the original emitted flux can 
be guessed from the signal analysis. 

As can be observed the analytical flux models do not fit perfectly the signal. Several 
fits seem actually rather poor (see for example TSNII-2). Nevertheless in almost all 
the cases in average the right value of the mass was obtained (see figs. I5.2H5.5|) . A 
bad fit to the flux model just increases the dispersion in the mass values, but does 
not change too much the best-fit point. Again, this is an indication that the right 
value of the neutrino mass can be measured even if the fit to the signal detailed 
shape is only approximate. 
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(a) TSNI-1 



(b) TSNI-2 





(c) TSNI-5 



(d) TSNI-6 





(e) TSNI-7 



(f) TSNI-8 



Figure 5.6: Fits to the flux model for various sets of input conditions using supernova 
model I. In each figure we compare the input time profile ft'^^{t) \5.^) (thick continuous 
line), the best- fit flux model ft^^^^^{t) for each signal (bunch of curves) and the observed 
time profile ft^^{t) (the outer and inner histograms enclose the observed profiles for all 
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(a) TSNII-1 



(b) TTSNII-2 
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(c) TSNII-5 



(d) TSNII-6 




(e) TSNII-7 



Figure 5.7: Fits to the flux model for various sets of input conditions using supernova 
model II. 



CHAPTER 5. RESULTS AND DISCUSSION 



5.3. QUANTIFYING THE SENSITIVITY OF THE METHOD 93 



The more flexible flux model II fits better the signals when supernova model II is 
used (compare TSNII-1 and TSNII-2). This flux model allows to describe more 
accurately the inflections of the flux for this model of supernova neutrino emission. 
The supernova model I signals are relatively well fitted by the simpler flux model I, 
and no significative difference in the best-fit mass can be observed with respect to 
flux model II. 

The parts of the detected flux that are fitted better are the initial rising phase and 
the decay. In all cases, when the model fits properly these phases, the fit to the 
mass is successful. 

Altogether, these evidences support the idea that the flux model fit is to a large 
extent uncorrelated with the neutrino mass measurement. Regardless of the avail- 
ability of a good astrophysical description of the signal, the method could extract 
enough information on the neutrino mass from a high statistics signal under a wide 
range of conditions. 

5.3 Quantifying the sensitivity of the method 

The results presented in the previous section give us a general idea of the overall 
properties of our method. Now we want to put quantitative limits on how much 
information on the neutrino mass could be extracted from a future supernova signal. 

In order to quantify the potential sensitivity of our method we have devised two 
different approaches: 

Averaging the upper bounds. If the neutrino mass is small and the sensitivity 
of the method is not large enough to resolve it, the only information that we can 
obtain is what is the largest neutrino mass compatible with the signal. In terms 
of the profile likelihood L the mass upperbound can be computed using the 
prescription ()4.10|) : 



We have performed a MC analysis generating ~ 40 synthetic signals per each set set 
of input conditions and assuming a MC neutrino mass mj^^ ~ 0. For each signal 
the mass upperbound with a 90% of probability was computed. We have used a flat 
prior probability for ^ and vanishing probability for ^ 0. 

Using these results we characterize the sensitivity of the method with two numbers: 
the average value of the upper bounds, and the dispersion Amup around this 




(5.3) 
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value (see fig. 15. 8|) . Clearly, the only meaning of these two numbers is that of 
identifying the possible range for the upper limits that could be obtained in a real 
case. In particular, Am^p is given only to indicate to what extent mup is a good 
representative of the different results of the entire ensemble, and should not be 
understood as the error on the quoted upper limit. 

Rejection of the massless neutrino hypothesis. If neutrino mass were large 
enough to be resolved, an alternative way to quantify the sensitivity of the method 
is to determine for which value of the input mass m^'^ the massless hypothesis can 
be rejected for a significative fraction of the analyzed signals. 

We performed different MC analysis using several values of m^'^ > 0. In each case 
for every signal we computed the lower limit using the prescription ()4.9p : 

/ p{ml\D,I) dmlc^ jC{ml\D,I) dml = CL (5.4) 

J m? J m? 

low low 

The minimum value of the input neutrino mass m^'^ = mmin for which m^^^ was 
larger than (90% CR) for more than 50% of the signals was used as a measure 
of the sensitivity (see fig. 15. 8p . In this case C must be normalized over the whole 
interval — oo < < oo to obtain a physically significant mi^^. 

In table l5.l1 we summarize the results of applying the previous procedures to quantify 
the sensitivity of the method. 

The first thing we can deduce from the results is that the average number of events 
in the signals is the dominant factor. This is clearly illustrated in fig. 15.91 were we 
have plotted the values of mup against the average number of events. It must be 
noticed that the improvement in the sensitivity due to an increase in the statistics 
is not dramatic. Comparing the sensitivity of the method at HK with the analogous 
results of SK we see that an increase by a factor of 10 in the statistics improves the 
sensitivity only by a factor of 2. 

Comparing the results of the combined analysis SK+KL (TSNI-5 and TSNII-5) 
with these performed on SK alone (TSNI-1 and TSNII-1) we see that sensitivity is 
determined by the larger statistics SK and the better energy resolution and lower 
threshold of KamLAND does not impact too much the result. A similar conclusion 
is obtained observing that the planned scintillator detector LENA presents only a 
minor improvement in the sensitivity with respect to SK. Even if the energy resolu- 
tion is better, the comparable fiducial volumes that implies comparable number of 
events results in quite similar sensitivities for the two detectors. 



CHAPTER 5. RESULTS AND DISCUSSION 



5.3. QUANTIFYING THE SENSITIVITY OF THE METHOD 95 



Upper bound intervals at 90% CL i 



::::j^p 



6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 
Sample Number 



o 




2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 
Sample Number 



Figure 5.8: Two ways to quantify the sensitivity of the method. Upper panel: the average 
of the mass upper bounds and its dispersion Afn^p. Lower panel: increasing the MC mass 
used to generate the samples we find the minimum value for which more than 50% of the 
analysis can reject the zero neutrino mass hypothesis. 
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1 10 100 

N. events ( x 10^ ) 

Figure 5.9: Relation between the number of events in the signal and the sensitivity of the 
method. The continuous line highlights the empirical "logarithmic" relation between the 
sensitivity and the amount of events in the signal. For different distances also the different 
TOF is involved in the determination of the sensitivity, and the respective results deviate 
from the linear-log fit. 
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Table 5.1: Results for rn^p and \J'rn^^^^ under different emission and detection conditions 
computed inside a 90% C.R. (95% C.R. values are in parenthesis) . The average number 
of events of the signals for each set of conditions is included for reference in the fourth 
column. The values reported in columns 5 and 6 are affected by statistical uncertainties at 
the level of ~ 5%. The typical number of signals per test is ~ 40. 



supernova model I 



Test 


Detector 


Distance 


N. events 


± Am^p (eV) 


V^iin (eV) 








(xlO^) 






TSNI-1 


SK 


10 kpc 


10.0 


1.0(1.1) ±0.2 


1.0(1.1) 


TSNI-3 




5 kpc 


40.0 


1.1(1.2) ±0.2 


1.1(1.2) 


TSNI-4 




15 kpc 


4.4 


1.3(1.4) ±0.3 


1.4(1.5) 


TSNI-5 


SK+KL 


10 kpc 


10.4 


1.0(1.1) ±0.2 


0.9(1.0) 


TSNI-6 


LENA 


10 kpc 


12.6 


0.9(1.0) ±0.2 


0.9(1.0) 


TSNI-7 


HK 


10 kpc 


170 


0.4(0.5) ±0.1 


0.5(0.6) 






supernova model II 




Test 


Inputs 


N. events 


m^p ± Am^p (eV) 


V^iin (eV) 








(xlO^) 






TSNII-1 


SK 


10 kpc 


5.9 


1.1(1.2) ±0.3 


1.1(1.2) 


TSNII-3 




5 kpc 


23.7 


1.1(1.2) ±0.3 


1.1(1.2) 


TSNII-4 




15 kpc 


2.6 


1.6(1.7) ±0.6 


1.6(1.8) 


TSNII-5 


SK+KL 


10 kpc 


6.1 


1.1(1.2) ±0.3 


1.1(1.2) 


TSNII-6 


LENA 


10 kpc 


7.5 


0.9(1.0) ±0.3 


1.1(1.2) 


TSNII-7 


HK 


10 kpc 


100 


0.5(0.6) ±0.1 


0.5(0.6) 


Reference test: 


10 kpc 


9.6 


0.8(0.9) ±0.2 


0.9(1.0) 
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Changing the supernova distance two competitive effects arise. On one hand de- 
creasing the distance increases the signal statistics, which in principle implies an 
increase of the sensitivity. However, at smaller distances the TOF delays induced 
by the mass are also smaller and become harder to identify. As can be learned with 
tests TSNI-3 and TSNII-3 these two effects tend to compensate each other. 

It is natural to ask if anything better could be done to measure neutrino masses 
from a supernova neutrino signal. In the attempt to answer this question, we have 
performed the following test: we have produced a sample of synthetic signals assum- 
ing no mixing in the spectrum (no oscillations) and using as inputs to our MC flux 
model I ()4.13|) with a suitable choice of the relevant parameters, together with an 
analytical a-distribution spectrum corresponding to the average energy profile given 
in fig. 15.1b . We have then performed our usual set of fits to the neutrino mass (we 
assume the SK detector) but fixing the value of the fiux parameters to the ones used 
in the MC, and we also adopt the same time varying spectrum used to generate the 
sample. This simulates the ideal (and unrealistic) situation where the time-energy 
dependence of the signal at the source is known, and the only free parameter is the 
neutrino mass. The results of this test are given in the last row in table 15.11 and 
should be compared with the results for tests TSNI-1 and TSNII-1. We see that only 
a mild improvement is achieved with respect to the realistic situation. This allows 
us to conclude that the sensitivity to neutrino masses of the detectors presently in 
operation is very likely bounded to values not much below 1 eV. 

5.4 Additional remarks on the method 

To arrive to the results we have just presented, some specific choices about the 
overall procedure that was followed had to be made. Here we will briefiy discuss the 
impact of some of these choices on the final results. 

Different priors. As was stated in sect. 14.4.21 the choice of the prior probability 
for mf, is a subtle issue. For all the analysis presented in this work we have used a 
function as the prior probability on m^. In order to understand the impact that 
a different prior could have on the results we have compared the mass limits with 
the ones obtained with a flat prior on m^. 

In flgure 15.101 we present the mass upper bounds obtained with both priors. We 
observe that the differences are not large. In the context of Bayesian reasoning this 
means that the signal contains enough "evidence" about the mass and therefore a 
change in the prior does not affect the results. 

Estimation of the energy spectrum. The choice of the spectral function used 
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Figure 5.10: Upper bounds obtained when two different prior probability on neutrino 
mass are used. Continuous line correspond a flat prior on ml as used for all the analysis 
reported in this work. Dashed line are the upperbounds for the same signals obtained with 
a flat prior on m^. 



to fit the energy distribution of tlie events represent an interesting issue in our 
analysis. In ref. [12] we used a Fermi-Dirac distribution. Comparing these results 
with the ones presented here, where the a-function was used, it is evident that the 
particular spectral shape parametrization has no big impact on the final results. 
The main difference is the more direct way the parameters of the a-distribution can 
be estimated. In the Fermi-Dirac case the estimation of the spectral temperature 
and effective degeneracy parameter rj involves a numerical procedure which mainly 
slows down the numerical analysis. 

Energy threshold. As was observed in sect. 15. 31 an increase in the energy resolution 
of the detector does not affect too much its sensitivity. However a too large energy 
threshold could affect considerably the method resolution. 

In ref. jHS] we compared the results for analysis performed on signals detected in SK 
assuming two different energy thresholds: 5 MeV, the threshold used here, and a 
more conservative threshold of 10 MeV. In that work the sensitivity dropped down 
by a factor of two making apparent that the low energy neutrinos play a central role 
in the determination of the sensitivity. 
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A new question arises: why the analysis of scintillator detector signals with a signif- 
icantly lower threshold does not produce a similar improvement in sensitivity? The 
answer is that the increase in number of events is about 1-2% (see fig. 13 .71) when 
we go from the 10 MeV threshold to 5 MeV (and this represents 100-200 additional 
low energy events in the SK signal) while lowering the threshold down to ~ 3MeV 
implies only a ~ 0.05 — 0.1% increase (5-10 events). 
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5.5 Summary and Conclusions 

We have described in this work a new method to measure and constrain the absolute 
scale of the neutrino mass using a high statistics signal from a future Galactic 
supernova. 

The method relies on three basic conditions: 

• An almost thermal neutrino spectra at the source (see sect. 12. 2() . 

• The distribution of the high energy neutrino events, less affected by the TOF 
delay induced by a non-zero neutrino mass, can be used to extrapolate the 
signal profiles at low energies, since the time scale of the spectral evolution is 
larger than the typical mass induced time delay (fig. 12. 8|) . 

• A theoretical description based on general characteristics of the neutrino flux 
time evolution (an early fast rise followed by a steady decay on time scales of 
several seconds) is sufficiently accurate to construct a likelihood for studying 
the neutrino mass more probable values. 

We combined a ML estimator and Bayesian techniques to construct the statistical 
procedure designed to constrain a neutrino mass, irrespectively of the particular flux 
model parameters used to describe the signal. 

Different tests corresponding to two different supernova models, different oscillation 
schemes and detector were carried out. For each particular set of conditions, about 
40 complete analysis were performed. 

The general results of these tests have been described in detail in sects. 15.21 In 
particular, it was shown that regardless of the fine details of the signal, the measured 
value of the neutrino mass can be nailed around the correct value. The spread in 
the uncertainty depends however, on the quality and the amount of details used to 
describe the signal. 

The two analytical models introduced to describe the flux behaved relatively well 
when applied to very different numerical neutrino fluxes. Satisfactory results where 
obtained with the first more simple flux model (|4.13p especially when fitting super- 
nova model I signals. The more flexible flux model II ()4.15p required a tuning of 
various parameter that had to be fixed to suitable value to leave just a reasonable 
number of free variables in the likelihood multiparameter extremization. 

The sensitivity of our method was estimated in two different ways, first, by deter- 
mining the typical upper bound on mj, that could be obtained in case the neutrino 
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mass is too small to be resolved. Secondly, by evaluating the minimum mass that 
could be distinguished from zero. 

We believe that the method that we have proposed represents an improvement with 
respect to previous techniques, both in sensitivity and in the independence from 
particular astrophysical assumptions. 

The sensitivity to neutrino masses of the detectors presently in operation can reach 
a level down to 1 eV. This is sizable better than present results from tritium /5-decay 
experiments fI7\ I28j. is competitive with the most conservative limits from neutri- 
noless double /3-decay [2312112111201, and although less precise than cosmological 
measurements [SBl OH IHE], is also remarkably less dependent from prior assump- 
tions. Future megaton water Cerencov detectors as Hyper-Kamiokande will allow 
for about a factor of two improvement in the sensitivity. However, they will not be 
competitive with the next generation of tritium /3-decay |94| and neutrinoless 
double /9-decay experiments (see and references therein). 

In figure l5?TT] we compare the limits that could be obtained with the present method, 
with other available and foreseen limits from laboratory, astrophysical and cosmo- 
logical data. 

We conclude that the occurrence of a Galactic supernova explosion within the next 
few years might still provide valuable informations on neutrino masses that however 
will not be able to explore a region much below 1 eV. Therefore as new laboratory 
experiments and cosmological observations will push the neutrino mass limits sen- 
sibly below 1 eV, the corresponding effects of the neutrino time of flight delays on a 
supernova signal will become unmeasurable. 
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Figure 5.11: Mass limits from present and future laboratory (MAINZ\^, KATBJNJJJ^. 
Oi^PPjMl), astrophysical fSNlM. EH \^) and cosmological studies (WMAPjMl)- The 
thick gray lines in the background indicates the absolute upper limit from Tritium decay 
rUi, < 2.3 eV \27l and the lower limit from oscillation evidences ruy > Am^^^^ ~ 0.05 
eV^. 
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Bayesian Inference: basic 
definitions 

The method and techniques used in this work to obtain informations about a neu- 
trino mass from a synthetic supernova signal are based on Bayesian principles of 
inference. Bayesian reasoning in data analysis has recently gained more and more 
relevance in a wide range of scientific disciplines, including frontier physics |97j . 
This appendix is aimed to introduce the main definitions and results of Bayesian 
inference in data analysis on which the method proposed in this work is based. A 
self contained and physics oriented introduction to Bayesian reasoning can be found 
in |3n] while a more complete review of Bayesian techniques and their applications 
in physics data analysis is given in pT7j . 

A.l Basic principles 

Bayesian inference is founded on two basic ideas. The first is the conception of 
probabihty as the degree of belief thai a given proposition is true as opposite to the 
conventional statistical definition of probability as the long-run relative frequency 
with which the event defined by the proposition occurred on many repeated exper- 
iments IHHl- Curiously (or indeed not) the Bayesian idea of probability correspond 
to the original concept by Laplace (as the father of inference reasoning in science) 
two centuries ago|98j. The second idea is that every probabilistic statement that 
we can make on any proposition must be a conditional probability, i.e. it must be 
conditioned to the available relevant information related with the proposition. 

Formally, Bayesian inference is based on two central rules that can be derived from 
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the requirement of logical consistency of the probability theory j^. By denoting 
with p{A\B) the probability that A is true, given that B is true, it is straightforward 
to write down the product rule: 



which states that the probability that propositions M and D are true given the 
background information / is equal to the probability of M, given D and /, times 
the probability of D given /. On the other hand we have the sum rule: 



where Mj is a set of mutually exclusive and exhaustive set of propositions and 1 refers 
to the convention to assign this value to the probability of a tautological proposition. 

A. 2 Bayes theorem 

Starting from the two rules introduced in the previous section two central results 
can be obtained: the Bayes theorem and the marginalization procedure. 

Let us suppose that in a given measurement we obtain the data D. Based on some 
background information / we assume that model M can describe the data. Using 
the previous knowledge (product of our expertise) on the phenomenon we assign to 
model M a prior probability p{M\I) which measure the plausibility of that model 
given the background information I. The posterior probability of M given the data 
and the background information can be computed from p{M\I) using the product 
rule ()A.ljl : 

v(D\M I) 

p{M\DJ)=p{M\lf-^-^^ (A.3) 

This result is called the Bayes theorem. Here p{D\M,I) is the probability that the 
data D be described by model M, and it is called the sampling probability for D 
or the Likelihood for model M. p{D\I) is called the evidence for D and represents 
the probability that the measurement produce the data D for the entire class of 
hypotheses (models). 

When M is described by a (continuous) set of parameters collectively denoted as 
{^}, posterior p{{6}\D,I) becomes multivariate probability distribution functions 



p{M,D\I) 



p{M\D,I) X p{D\I) 
p{D\MJ) X p{M\I). 



(Al) 




(A.2) 
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(pdf) of the parameters, while the Likehhood p{D\{6},I), that we will denote by 
the symbol C{D; {6}) in spite of its explicit dependence is not by itself a pdf for the 
parameters. The evidence p{D\I) is independent of {6} and plays simply the role of 
the posterior pdf normalization constant = p{D\I). Bayes theorem now reads: 



where for simplicity we have dropped out the reference to the background informa- 
tion I. 

An interesting way to see what Bayes theorem establishes is to think that our knowl- 
edge on how a model describes a given phenomenon grows when new evidence is 
accumulated. In this sense we could thing at Bayes theorem as a recipe for learning 

nnoi. 



A. 3 The marginalization procedure 

Often some of the parameters that are used to describe a model, thought important 
for the computation of the Likelihood, are uninteresting for the final conclusions that 
one wants to extract from the experimental evidence. The principles of Bayesian 
reasoning provide a natural way to marginalize these nuisance parameters. Using 
the continuous limit of the sum rule (IA.2I): 



jp{9\D)d9 = l, (A.5) 
the marginal posterior probability for parameter 6*0 can be written as: 



p{eo\D) = j d{e}p{eo,mD) 

= N~^p{eo) j d{e}c{D-e,,{e})pm) (a.g) 

Marginalization is one of the most important features of Bayesian inference when 
compared with conventional approaches. However, performing the multidimensional 
integrals required for finding the marginal distribution can be a computational chal- 
lenge in terms of the large times necessary to carry out such integrations. Several 
solutions have been devised and are used in the most complicated problems where 
Bayesian inference is applied (see section 2 of ref. |lUUj ). As was described in sect. 14. 6( 
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we have used the simple approximation of taking the profile likelihood (PL) as the 
marginal distribution, under the assumption that the likelihood can be approximated 
by a multivariate Gaussian (normal) pdf. In the following section we will prove that 
the profile of a normal likelihood coincides with its marginal distribution. 

A. 3.1 Marginalization of a normal pdf 

The general form of a multivariate normal pdf for parameters 9^ = {9Q...9n~i)'^ is: 



^"^^"^^^'^^= (2.)-/Vdet(C) ^"^ 



.i(r-/r)Tc-i(r-/r 



(A.7) 



where /i = {^o...fin-i)^ are the mean values of the parameters (here for simplicity 
we will assume /ij = 0) and C is the covariance matrix: 

O'o PoiCToO"! •• 
C = I poi(To(Ti a\ ... I (A.8) 

with pij the correlation coefficient between parameters i and j. 

To simplify let's take the particular case of a bivariate normal pdf. In this situation 
eq. ()A.7|1 is given by: 



^2(6'o,6'i;cro,a-i) 



1 



exp 



exp 



2(1 -P^ 
[ael^hel - 2c Ml) 



2 ' 2 



(ToCTi 



(A.9) 



with p = pi2 and in the last equation we have defined Z 
1/aiil - p2), b = (1 - p2) and c = p/aoa^{l - p'). 

Using ()A.9j) the marginal pdf on 6*0, p{0o) reads, 



p{0o 



d9iG2{0o,9i;(To,(Ti) 
^ V 2 26 



dOi exp 



'-M9^ + df 
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exp 



2(T, 



2"0 



27rcro 



exp 



2ai 



where after completing the square in the exponential of the second equation [d = 
c9o/Vb) the resulting integration is performed using J exp(— x^/2(t) = \/27ra. As 
can be see from the last equation, the marginal distribution in is simply a normal 
distribution with dispersion ao- 

Now we will proceed to compute the profile likelihood (PL) p{6o) which is defined 
by: 



p{9o) OC g{9o, 6'imax; O-Q, CTl) (A. 10) 

where ^imax refers to the value of 6i for which Q is maximum at the given value of 
^0 (see figure WTT^ . The proportional symbol 'oc' means that p{6o) is not necessarily 
normalized. 

The value of 6'imax is found maximizing the denominator in the exponential of ()A.9j) 
which gives 6'imax = cOq/B. This results in, 



^(6^0) ~ exp 
~ exp 



2 bc^ 2 2 
2 v°^o + -^^o-2y^o' 



'2cr2 ° 



(A.ll) 
(A.12) 



Therefore the profile likelihood is also a normal distribution in ^0 with dispersion 
equal to ctq and coincides with the marginal distribution p{6o)- 

This result can be generalized to more than two variables. When the likelihood is a 
multivariate normal distribution in the parameters the profile likelihood for a given 
parameter coincides with the marginal distribution. 



A. 4 Credible regions 

(Marginal) posterior pdf are the general outcome of inference analysis based on 
Bayesian reasoning. It is common to summarize the properties of these distributions 
using for example the mode (the values of parameters where the pdf is maximum) 
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or the mean value and dispersion of the parameters as obtained from the posterior 
using the classical prescriptions. 

It is also common to construct credible regions (CR), defined as the interval (region) 
where the parameter value (s) is found with a probability CR, 

I dep{e\D) = CR (A.13) 

JcR 

The credible regions are the analog of the confidence intervals in frequentist statistics 
that however are computed and interpreted in a very different way. 



APPENDIX A. BAYESIAN INFERENCE: BASIC DEFINITIONS 



Appendix 



SUNG: supernova Neutrino 
Generation tool 

In chapter El we describe how the properties of a supernova signal are computed 
starting with the characteristics of the emitted neutrino signal as provided by self- 
consistent simulation of the core collapse and using the proper description of neu- 
trinos oscillations in the supernova mantle and the Earth interior and simulating 
the neutrino detection process. Numerically this procedure involves a particular 
manipulation of the supernova simulation results that are used as the inputs of the 
general description of the signal and the development of specially designed routines 
to model neutrino oscillations under general circumstances in the supernova mantle 
and the Earth interior. 

Other numerical challenges arise when the purpose is to generate full statistics sig- 
nals. In this case it is necessary on one hand to use especialized algorithms and 
techniques to sample the complex observed rate but in the other hand to design 
mechanisms to perform a cross-checking to the generated signal to ensure that its 
events already follows the original distribution. 

Using the experience and the computer codes created in the frame of this project 
we have created a computer package, SUNG (SUpernova Neutrino Generation tool) 
intended to offer a simple and acsequible computer solution to reproduce the results 
published in this work and to perform similar aproaches to the analysis of supernova 
neutrino signals. 

For new or experiencied researchers in the field SUNG is a solution to fastly obtain 
results that can require non-trivial programming efforts. Certainly it is not an ul- 
timate solution for the computation and generation of supernova neutrino signals 
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under rather arbitrary conditions. The main philosophy behind this effort is that 
of promote the creation of computational solutions in the field of neutrino astro- 
physics as has been done in other fields (e.g. CMBFAST for CMBR simulation and 
analysis |lUlj ). 

In this appendix we describe the general features of SUNG and some of the methods 
and algorithms used to perform the most important tasks that the package can 
perform. 

For a more complete description of the package including simple reference manuals 



for the final user and the programmer please refer to the web site http : //www . sungweb . tk 



B.l Package structure 

SUNG is more than a simple single program to generate supernova neutrino signals. 
It is also a especialized programming framework that can be used for a wide variety 
of different applications. 

SUNG has three main components: 

• A library of (ANSI C) routines for data analysis originally designed or adapted 
from other libraries, and other routines created to simulate neutrino oscilla- 
tions in the Supernova mantle and inside the Earth. 

The routines has been designed and organized to offer a coherent programming 
framework. Even the especialized numerical routines of other scientific libraries 
are inserted into routines especially formated for this library jlU2j . 

In the core of the numerical integration, special functions, interpolation and 
random number generation routines SUNG uses routines of the GNU scientific 
library fGSL) [Tn5] . Function minimization routines used in fit procedures use 
the MINUIT package routines from the CERN hbrary [TUlj . 

The SUNG library also includes a series of routines to generate data plots in 
eps and ps format using the plotting package GNUPLOT |lU5j . 

In this sense the SUNG library can be seen as a special purpose common 
interface to those scientific libraries and packages. 

• A set of C programs, written down using the routines in the SUNG library and 
perl/tcl scripts designed to perform mainly four tasks: manipulate super- 
nova simulation results, compute the properties of supernova neutrino signals 
under different input conditions, generate full-statistics realizations of those 
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signals and analyze high-statistics supernova neutrino signals using a suitable 
analytical model. 

All the programs are designed to be ran from the command-line. That pro- 
grams are the best example on how complex tasks can be programmed using 
the routines of the package library. 

The package have a special Graphical User Interface (GUI) which has been 
designed to run on a web browser. This webGUI consist of an elaborated set of 
php scripts and works almost as a regular off-line GUI but with the advantage 
of not depend on any particular GUI library. Once installed in a server (with 
an Apache and Php server) it can be run locally or remotely on a platform 
independent basis. 

A running version of the SUNG web interface can be found in the web site of 
the project http://www . sungweb . tk 



B.2 Methods and Algorithms 

In the following paragraphs we describe the methods and algorithms used by the 
package to perform the most important general tasks for which it was designed. 

It is important to recall that this work was performed using the same computer 
codes of the package and therefore the algorithms described below are the same 
ones involved in the production of the results presented here. 

B.2.1 Manipulation of Supernova simulations 

In order to describe the properties of an observed supernova neutrino signal it is 
necessary to know the properties of the emitted neutrinos. 

There are three main pieces of information that describe almost completely the 
neutrino emission from a supernova at a given time: the luminosity, the average 
energy and the non-thermal distorsion (quantified with an a-parameter for example). 
In the ideal case the non-thermal distortions information is replaced by the detailed 
knowledge of the spectral shape at each time. 

Most of the supernova simulations results produce as an output that informations 
and in our case they are the input of the computation of the observed signal prop- 
erties. 

In the manipulations of results from supernova simulations several problems must 
be faced. One of the most common problems arise when a given simulation result 
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is limited to a shorter interval than the total assumed duration of the signal (20 sec 
in this work). This is the case for Supernova model II [BJ where neutrino emission 
properties are just reliable for times less than 300 msec. 

The routines included in the package are not prepared to face this kind of limitations 
and an external extrapolation of the results is required (see sec. 15. ip . 

All the data files describing the evolution and the luminosity of the spectral prop- 
erties are treated by the SUNG library as continuous functions using linear interpo- 
lation. 

B.2. 2 Expected signal properties 

The signal properties computed by SUNG are the same described in sec. 13.31 
• Number of events. This property is computed using the definition ()3.46|) : 

N{AE^,,Atu) = PdE ^dt^^^^^. (B.l) 
Je, Jt^ dE dt 

Where the detected signal rate is given by ()3.4H) and ()3.40|) : 



d^nT,^{E,t) 
dEdt 

L^St{E,t) = pS^+{l-p)S^. 



Nt j dE' S^'\E', t) a{E') e{E') n{E, E') 



The three dimensional integration involved in the computation of is per- 
formed using a recursive procedure |106j supported by one- dimensional adap- 
tative Gaussian-Kronrod integrations. 

Time-profiles. There are two kind of time-profiles computed by SUNG. The 
energy integrated time-profile as defined in sec. 13.3.31 also called in the package 
context the marginal time-profile (j3.49j) : 

/r(t)= fdE^^^. (B.2) 
^ ^' J dEdt ^ ' 

The other profile is given by the histogram that describes the number of events 
inside narrow time windows of regular width At. This histogram is given by: 
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ti-At/2 



dE dt 



Energy-profiles. The energy-profiles (time integrated spectrum and energy- 
profile histogram) are computed similarly using the definitions given in sec. 13.3.1^ 



B.2.3 Signal generation 

The generation of full statistics signals is the central problem solved by SUNG. 

As explained in sec. 13.41 we can conceive several alternative algorithms to generate 
the events in a supernova neutrino signal including all the effects of the neutrino 
propagation, oscillation and detection. 

SUNG uses the approach of generate the signal by sampling directly the detected 
rate. 

The algorithm used by SUNG to perform this task is described in the following lines: 



Signal generation Algorithm 



1. Compute the number of events that different neutrinos emitted as 
a given flavor will produce in the detector: 



-1 'J^e^l^e 

where N-p^_,-p^ i^Vx^uJ ^^e number of I/g events produced by 
neutrinos emitted from the supernova core as Vf. {Vr^). 

2. Reshuffle the total number of neutrinos computed in the previ- 
ous step according to a Poisson distribution with average equal 
to A^'tot = N-p^_^j!^ + N-p^^v^- Acoording to the reshuffled value A^'tot 
compute the finally detected number of neutrinos of each flavor 
-A^e,i7a:^Fe obcyiug the original proportions. 

3. Generate N^^^^^ (^NL^^-^) pairs of neutrino energies and times 
{Ei,ti\ according to the respective neutrino detected flux (|3.38p . 
pSe/L'' ((1 -p) S^/L^) with S-e,^ given by (USED: 



dE / dtNxp 



dE I dtNril -p)^ 
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Sa{E,t) 



Lait) 



(B.4) 



4. Shift the neutrino time according to the mass induced delay (jOJ: 



5. Reset the neutrino times in order to place the first neutrino detected 
at ti = 0. 

6. Compute the enegry of the secondary produced by the detection of 
each neutrino: Ef = Ei — Qpn- 

7. Reshuffle the positron energy according to the energy measurement 
uncertainty ()3.45p : 



The result of this procedure is the set {E'f, ti} with z = 1, . . . , A^'tot- 



ti=ti + Atto{{ml,L,Ei) 



(B.5) 



E'^^gix = Ef,a = AE) 



(B.6) 
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C.l List of symbols 



Symbol 


Definition 




Absolute scale of neutrino masses. In the case of non-degenerate 
masses it refers to the electron neutrino mass in the context of this 
work. 




Non-electron neutrinos and antineutrinos = t'^,r) ^^/^.t- 


L 


Supernova distance. 




Radius and mass of a neutron star. 


Attof 


Mass induced time-of-fiight delay. 


dTf^/dt 


Rate of neutrino emission, i.e. number of neutrinos of any energy 
emitted per unit of time. 


F{E) 


Neutrino spectrum at the source. 


P 


Pinching parameter, p oc (E"^) / (E). 




Neutrino flavor eigenstates. 
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Symbol Definition 



^Af(t^) Neutrino mass eigenstates in vacuum (in matter). 



Uni entries of the MNS matrix. 



In matter effective neutrino masses. 



7 Adiabaticity parameter. 

Ph,l {Ph,l) Supernova mantle jumping probabilities. 

Pai [Pai) Supemova mantle conversion probabilities = P(z/q, z/j). 



P^p (P-^-^) Conversion probabilities = Pji'a yp). 
P® (P®) Conversion probabilities in the Earth. 



So,{E, t) Total flux of neutrino flavor a at the source, i.e. number of neutri- 
nos a emitted per unit of time and per unit of energy. 



S^^^{E, t) Total flux of neutrino flavor a arriving to the detector, i.e. number 
of neutrinos a arriving per unit of time and per unit of energy. 



Ea Average energy of neutrinos a at the source. 

d^Ua/ dE dt Total rate of neutrino flavor a in the detector, i.e. number of neu- 
trinos a detected per unit of time and per unit of energy. 



fsiE) Time-integrated energy spectrum. 



ft{t) Energy-integrated time profile. 



Ai A supernova emission model. If the model can be parametrized 

with a set of parameters {6} the model is represented by A^(^). 
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Symbol 


Definition 


V 


Data set, i.e. set of energy and time pairs \Ei,ti\. 


f{E,t) 


Probability distribution function used to describe the signal (model 
density probability). 


r 


ijiKeiinooci iuncLion. 


m 


Flux model. 




) Normal distribution with mean /i and standard deviation a . 


f{E, t) 


Regularized model density probability. 


t 


Profile Likelihood. 
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C.2 List of abbreviations 



I\ Drt3 V . 


t3 n 11 1 L 1 o 11 


SN 


bupernova/Supernovae. 


TOF 


lime-oi-nignt. 


LTE 


Local i nermodynamic Lquilibrium. 


MNS 


Maki-iNlakagawa-Sakata matrix. 


NH 


AT 1 1 * 1 

iNormal hierarchy. 


IH 


T J 1 1 * 1 

inverted hierarchy. 


LF 


T ■ 1 1*1 1 r J • 

Likehhood lunction. 


log-LF 


T 'jI fT'll'l 1 T — ' J • 

Logarithm oi Likehhood Inunction. 


ML 


IV/T * T*11*1 1 1 

Maximum Likehhood analysis. 


A/TP 


ivioiiLe v^aiio. 


CR 


Credible regions (credible region probability). 


PL 


Profile likelihood. 


SK 


Super-Kamiokande. 


KL 


KamLAND. 


HK 


Hyper-Kamiokande. 
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